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ABSTRACT 


A  Collision  Tolerant  Pile  Structure  (CTPS)  for  deploying  navigation  aids  in 
shallow  water  was  developed  having  the  ability  to  sustain  collisions  by  barge 
traffic.  The  CTPS  concept  considered  consists  of  the  aid  itself  mounted  at  the 
top  of  a  rigid  pile  which  is  hinged  just  above  the  mudline.  The  hinge  is  omni¬ 
directional  and  allows  the  pile  to  fold  down  in  the  event  of  a  collision.  The 
hinge  also  provides  a  restoring  moment  to  return  the  structure  to  the  vertical 
position. 

The  CTPS  was  designed  to  meet  Coast  Guard  criteria.  The  pile  inclination 
angle  must  be  small  (5  deg  as  a  design  goal,  10  deg  maximum)  during  extreme 
operating  conditions  (specified  as  30  ft  depth,  3  kts  current,  60  kts  wind  and 
5  ft  waves).  The  system  must  survive  hurricane  conditions  (that  is,  9  ft  of 
storm  surge,  100  kts  winds  and  6  ft  waves).  The  CTPS  is  required  to  sustain 
impact  by  a  barge  moving  at  10  kts,  and  it  is  desired  that  the  mudline  clearance 
be  as  low  as  3  ft.  The  CTPS  installation  must  be  similar  to  existing  practice, 
and  cost  must  be  low  ($5,000  as  a  design  goal,  $10,000  maximum). 

To  serve  as  a  design  tool,  computer  programs  were  developed  for  simulating 
pile  dynamics  under  the  conditions  referred  to  in  the  design  specifications. 
Separate  programs  were  written  for  operating,  hurricane  and  collision  conditions 
and  for  the  recovery  process.  The  CTPS  was  modeled  as  a  mass- rigid  bar- hinge 
system  with  loadings  and  resulting  motion  confined  to  a  vertical  plane.  Wind 
is  considered  steady;  wave  activity  is  represented  by  a  regular,  sinusoidal 
wave,  and  current  is  taken  to  be  steady  and  uniform  with  depth.  During  the 
collision  process  the  barge  is  assumed  to  maintain  constant  speed.  In  use, 
design  and  environmental  parameters  are  input,  and  the  program  calculates  the 
response  time  series. 


Using  preliminary  computer  model  results  for  guidance,  the  hinge  component 
was  developed.  Several  hinge  concepts  were  evaluated,  and  a  central  universal 
joint,  peripheral  stay  arrangement  was  determined  to  be  the  only  satisfactory 
system.  A  major  feature  of  this  concept  is  the  use  of  a  single,  pre-stressed 
spring,  housed  within  the  pile  (made  of  pipe),  to  tension  the  stays.  A  design 
based  on  this  concept  was  developed,  and  a  physical  scale  (1/15)  model  was  built 
for  testing. 

Initial  experiments  were  conducted  out  of  water  to  determine  the  hinge 
stiffness  characteristics.  These  preliminary  results  indicated  that  the  stiffness 
is  adequate  for  meeting  the  vertical ity  requirement  and  for  recovery  from  a  knock¬ 
down.  Next,  tests  were  conducted  at  the  U.S.  Coast  Guard  Academy's  circulating 
water  channel.  Here  it  was  found  that  the  scale  model  CTPS  possessed  more  than 
sufficient  stiffness  to  meet  the  verticality  requirement  at  (Froude  scaled) 
maximum  current.  Hydrodynamic  load  measurements  showed  that  the  drag  coefficient 
was  less  than  one  thus  justifying  the  use  of  unity  in  the  computer  simulations. 
Collision  experiments  were  then  carried  out  in  the  UNH  indoor  pool.  It  was 
found  that  the  scale  model  could  endure  typical  barge  collisions  at  the  design 
speed  (Froude  scaled)  and  return  to  the  vertical.  These  observational  results 
were  documented  by  taking  movies  and  photographs.  Impact  force  measurements 
confirmed  that  the  computer  simulations  were  sufficiently  accurate  (within  23%) 
for  design  purposes. 

To  estimate  full-scale  CTPS  dynamic  response,  the  computer  models  were 
applied  to  a  full  scale  prototype  under  the  design  criteria  conditions.  It  is 
predicted  that  the  maximum  inclination  angle  restriction  (10  deg)  can  be  met  and 
that  no  damage  to  the  pile  and  hinge  should  occur  during  hurricane  conditions. 

The  pile  system  should  survive  typical  barge  (12  ft  draft)  collisions  at  maximum 
speed,  but  impact  loads  were  found  to  increase  enormously  with  decrease  in 
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clearance.  For  this  reason  and  because  height  is  needed  for  the  hinge  itself, 
the  3  ft  mudline  clearance  specification  cannot  be  met. 

Logistical  considerations  were  investigated,  and  procedures  for  installing 
the  pile  using  familiar  methods  were  outlined.  Cost  estimates  indicate  that  the 
CTPS  itself  meets  the  maximum  allowable  restriction  ($10,000)  and  that  CTPS 
systems  should  become  cost  effective  after  3  collisions. 

Because  nearly  all  the  design  specifications  have  been  shown  to  be  achievable, 
it  is  concluded  that  the  CTPS  design  should  be  successful.  We  therefore  re¬ 
commend  that  the  next  step  in  development  be  taken  in  which  a  full  scale  proto¬ 
type  is  built  and  field  tested. 
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I.  INTRODUCTION 


BACKGROUND 

Rigid  pile  structures  currently  used  by  the  Coast  Guard  to  support  aid-to- 
navigation  markers  are  susceptible  to  collision  by  towed  barges.  Tugs  towing 
barges  skirt  the  edges  of  narrow  channels  to  avoid  deep-draft  ocean  vessels 
operating  in  the  main  channel.  Though  the  tugs  themselves  can  usually  avoid 
navigation  aids  on  channel  boundaries,  the  barge  strings  being  towed  will  often 
hit  the  navigation  aid  supporting  structure.  Considerable  expense  is  involved 
in  locating  the  lost  position,  pulling  the  damaged  marker  and  driving  a  new  pile, 
as  well  as  purchasing  structure  components. 

This  problem  can  possibly  be  alleviated  by  replacing  existing  rigid  pile 
supports  by  a  compliant  Collison  Tolerant  Pile  Structure  (CTPS).  Miller  (1982) 
provides  evidence  that  this  concept  is  promising  and  should  be  developed.  This 
report  presents  results  of  a  study  to  develop,  analyze  and  test  a  CTPS  design 
able  to  sustain  multiple  collision  without  requiring  major  repairs. 

The  CTPS  concept  considered  consists  of  an  aid  to  navigation  marker/light 
mounted  on  a  single  pile  which  is  hinged  just  above  the  mudline  as  shown  in 
Fig.  1.  The  hinge  is  omnidirectional  with  respect  to  barge  course  and  should 
possess  full  vertical  to  horizontal  articulation.  The  hinge  must  provide  a 
restoring  moment  to  return  the  pile  to  the  upright  position  after  a  knockdown 
and  to  maintain  a  near  vertical  position  while  the  system  is  in  normal  operation. 
The  CTPS  system  and  the  hinge  component  in  particular  were  designed  to  meet 
specific  performance  requirements  formulated  by  the  Coast  Guard. 

DESIGN  REQUIREMENTS 

The  CTPS  Is  to  be  suitable  for  use  in  the  port  of  Houston,  TX,  an  area  for 
which  pile  destruction  is  especially  serious.  Design  requirements  include  the 


capability  of  being  installed  by  existing  Coast  Guard  construction  vessels 
which  can  routinely  drive  piling  but  do  not  have  crews  capable  of  sophisticated 
welding  or  metal  work.  Logistical  design  requirements  are  given  in  Table  1. 
Operational  requirements  include  maintaining  a  position  within  an  "envelope11  of 
near-vertical  inclination  angles  when  acted  upon  by  weight  and  environmental 
loads  such  as  wind,  wave  and  current  forces.  Design  verticality  limits,  weight 
and  environmental  loadings  are  specified  in  Table  2.  As  a  design  goal,  the 
structure  should  also  survive  (without  daymark)  the  hurricane  conditions  given 
in  Table  3.  Lastly  the  Coast  Guard  requests  that  the  CTPS  be  able  to  withstand 
the  barge  collision  outlined  in  Table  4. 

The  CTPS  design  life  must  be  5  years.  The  design  goal  for  cost  is  under 
$5,000,  while  the  maximum  cost,  based  on  10  units/year,  must  not  exceed  $10,000. 

DESIGN  APPROACH 

Determining  the  feasibility  of  meeting  the  design  requirements  and  developing 
the  best  design  solution  for  the  hinge  component  were  accomplished  using  computer 
simulations  and  laboratory  testing  of  a  physical  scale  model.  Major  elements 
of  the  design  and  analysis  program  include  computer  modeling  of  the  pile  dynamics, 
hinge  component  development,  construction  of  a  physical  scale  model,  testing 
of  the  model's  verticality  performance  in  a  water  channel  and  observing  the 
model's  collision  response  during  experiments  conducted  in  a  pool. 

The  initial  task  was  to  develop  computer  programs  to  predict  pile  performance 
under  the  various  conditions  referred  to  in  the  design  requirements.  Assumptions 
and  key  equations  in  the  mathematical  development  are  outlined  in  Section  II. 

The  computer  simulations  were  used  to  relate  the  overall  design  criteria  to  hinge 
stiffness  characteristics  thus  providing  a  basis  for  hinge  component  development. 
The  computer  models  were  also  used  as  an  aid  in  interpreting  experimental  data 
and  to  predict  the  performance  of  a  full-scale  prototype. 
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Table  1.  Logistical  requirements 


(a)  Installation,  maintenance  and  removal  of  the  structure  must  be 
accomplished  by  existing  Coast  Guard  construction  vessels  and 
must  not  require  the  use  of  divers  or  sophisticated  underwater 
equipment. 

(b)  Proposed  structure  must  not  impose  serious  departures  from  present 
practices  of  pile  fabrication  and  installation  and  must  fit  within 
the  capabilities  of  existing  vessels. 

(c)  Any  single  component  of  the  structure  must  not  exceed  50  feet  by 
20  inches  in  size  or  weight  more  than  10,000  lbs. 

(d)  When  removed  for  repair  or  replacement  the  designed  structure  must 
be  retrieved  intact  or  constructed  to  separate  at  or  below  the 
mudline  in  order  not  to  become  a  hazard  to  navigation. 


Table  2.  Vertical ity,  weight  load  and  environmental  requirements. 
The  structure  must  support  the  dead  load  plus  live  load  and  two 
daymarks.  The  structure  must  meet  the  vertical ity  requirement 
with  dead  load  and  daymarks  under  the  environmental  conditions. 


Vertical ity 

+  5°  (design  goal) 

+  10°  (maximum  allowable) 

Dead  Load 

350  lbs 

Live  Load 

200  lbs 

Daymarks 

2 

2  at  36  ft  each  located  7 

ft  above  MHW 

Environmental  conditions  Design  goal 


Minimum  allowable 


Water  Depth  Maximum 
Minimum 


30  ft. 
10  ft. 


20  ft. 
15  ft. 


Current 


3  kts. 


2  kts. 


Wind 


Waves  -  Height 
Period 


60  kts  (gusts) 

50  kts  (sustained) 

5  ft. 

3-5  sec. 

15° 


60  kts  (gusts) 

40  kts  (sustained) 

4  ft. 

3-4  sec. 

15° 


Bottom  Slope 
Bottom  Consistency 


Soft  Clay 


Soft  Clay 


Table  3.  Hurricane  conditions.  As  a  design  goal  the  structure 
should  survive  and  return  to  normal  operation  after 
the  following  hurricane  conditions  with  dead  load 
less  daymark. 


Hurricane  Conditions 


Water  Depth  Maximum 
Minimum 


Storm  Surge 

Current 

Wind 

Waves  -  Height 
Peri od 


Quantitative  Extent 

30  ft.  (design  goal),  20  ft.  (allowable) 
10  ft.  (design  goal),  15  ft.  (allowable) 

9  ft. 

3  kts.  (design  goal),  2  kts.  (allowable) 

100  kts.  (gusts,  75  kts.  (sustained) 

6  ft. 

4-6  sec. 

15° 


Bottom  Slope 
Bottom  Consistency 


Soft  Clay 


Table  4.  Barge  Collisions.  The  structure  must  survive  the 
following  collision  conditions  and  return  to  normal 

operation. 


Collision  Condition 
Barge  speed 

Bottom  clearance  between 
mudline  and  barge 

Number  of  collisions 


Quantitative  Extent 
10  kts.  (maximum) 

3  ft.  (minimum) 

5  per  year  (maximum) 


Using  guidelines  established  by  the  computer  modeling  effort,  a  hinge 
concept  was  chosen  and  developed  as  discussed  in  Section  III.  The  selected  design 
and  construction  details  for  scale  model  fabrication  are  presented  in  Section  IV. 
Next  the  actual  hinge  stiffness  characteristics  were  measured  in  out-of-water 
experiments  (referred  to  in  this  report  as  "bench"  tests).  Results  for  hinge 
moment  as  a  function  of  angle  are  given  in  Section  V. 

The  water  channel,  vertical ity  tests  conducted  at  the  Coast  Guard  Academy 
and  collision  tests  carried  out  at  the  UNH  indoor  pool  are  discussed  in  Sections 
VI  and  VII,  respectively.  These  experimental  results  were  then  used  to  refine 
parameters  selected  for  the  final  set  of  computer  simulations,  described  in 
Section  VIII,  which  predict  full-scale  behavior  of  the  selected  design.  Logistical 
and  cost  factors  are  presented  in  Section  IX.  Conclusions  drawn  from  this  study 
are  summarized  and  recommendations  for  prototype  development  are  made  in 
Section  X. 


II.  PILE  DYNAMICS  MODELING 


MODELING  APPROACH 

Computer  programs  were  developed  for  modeling  the  dynamics  of  pile  systems 
under  conditions  referred  to  in  the  design  criteria.  Since  the  overall  CTPS 
design  criteria  given  in  Tables  2-4  refer  to  distinct  and  very  different  conditions, 
several  computer  models  were  developed.  Programs  were  written  for  operating 
conditions  in  which  small  angle  verticality  restrictions  must  be  met,  hurricane 
conditions  in  which  forcing  and  angular  motion  may  be  larger,  and  collision 
conditions  in  which  barge  contact  is  the  dominant  feature.  All  models,  however, 
share  some  common  assumptions  and,  as  a  consequence,  many  dynamic  equations  are 
the  same  for  all  applications.  General  features  of  the  modeling  approach  are 
discussed  here,  while  the  specifics  of  individual  computer  programs  are  detailed 
in  the  following  subsections.  Program  listings  are  given  in  Appendix  A. 

The  CTPS  is  considered  to  be  a  flexible  hinge-rigid  beam-mass  system  such 
as  that  shown  in  Fig.  2.  The  hinge  is  omnidirectional  and  possesses  restoring 
moment  stiffness.  Weight  and  current  forcing  are  external  loads  common  to  all 
major  computer  models,  while  wind,  wave  and  barge  contact  forcing  may  or  may  not 
be  present  depending  on  the  application.  In  all  models,  the  directions  of  current 
and  (when  present)  wind,  wave  and  barge  motion  are  assumed  col  linear  corresponding 
to  the  worst  case  situation. 

The  governing  dynamic  equation  for  the  hinge-beam-mass  system  considered  is 
the  time  rate  of  change  of  angular  momentum  equation  applied  at  the  (fixed 
point)  hinge, 

IH  0  =  I  M's  (1) 

where  IH  =  moment  of  inertia  about  the  hinge,  8  =  angle  of  pile  with  respect  to 
the  vertical,  (*')  indicates  two  derivatives  of  (  )  with  respect  to  time  t,  and 
M  refers  to  moments  applied  about  the  hinge.  (All  terminology  used  is  summarized 
in  the  NOMENCLATURE  section). 


CTPS  nomenclature.  The  CTPS  consists  of  a  flexible 
hinge,  a  rigid  beam  pile,  a  mass  near  the  top  representing 
the  load,  and  the  dajiarlc  boards.  Barge  and  environmental 
parameters  are  also  Indicated. 


The  flexible  hinge  will  be  constructed  to  provide  a  restoring  moment  to 
the  CTPS.  It  is  desired  to  have  the  hinge  be  very  stiff  at  small  angles  to  meet 
the  vertical ity  requirement  under  operating  conditions.  Yet  hinge  moment  stiffness 
at  the  large  angles  encountered  during  collisions  should  be  limited  in  order  to 
reduce  maximum  pile  bending  moment.  The  piece-wise  linear  behavior  of  hinge 
moment  -  MH(0)  shown  in  Fig.  3  has  these  characteristics  and  is  used  in 
the  models.  The  hinge's  behavior  at  small  angles  is  determined  by  the  initial 
stiffness  (slope)  k-j,  while  properties  at  large  angles  are  additionally  influenced 
by  the  large  angle  stiffness  (slope)  k2.  It  is  desirable  to  have  k-j  »  k2,  and 
the  breakpoint  angle  (point  of  slope  change)  should  be  at  the  limit  of  the  CTPS's 
operating  range. 

The  upsetting  gravitational  moment,  Mq,  due  to  pile  and  load  weight  has  the 
mathematical  form 

Mg  =  im  Sine  W,  *  1/2  !p  sine  Wp  (2) 

where  lengths  2^  and  2p  are  shown  in  Fig.  2  and  W £  and  Wp  are  load  and  pile 
weights,  respectively.  Hollow  piles  are  assumed  free- flooding,  and  the  restoring 
moment  effect  due  to  bouyancy  is  neglected. 

The  moment  load  induced  by  relative  water  movement,  M^,  is  evaluated  using  a 
form  of  Morrison's  equation, 

As  d  2 

MC  =  fQ  sCI  pw  Cw  dpur2  +  Cm(JL4_)  pw  Gr]  ds  (3) 

where  s  =  pile  coordinate  shown  in  Fig.  2,  =  submerged  length,  pw  =  water 

density,  Cw  =  drag  coefficient  of  the  pile  in  water,  =  inertia  coefficient  of 
the  pile,  dp  *  pile  diameter,  and  up  =  relative  velocity  normal  to  the  pile.  The 
relative  velocity  includes  steady  current,  motion  of  the  pile  itself  and  wave  fluid 
velocity.  Wave  motion  is  taken  as  that  of  a  regular  (single  frequency),  small 
amplitude  (linear)  surface  wave.  Thus  the  wave  fluid  velocity  components  are 


)  cos(kx-ot) 


oHw  cosh  k(dt+y) 
uw  T  ^  sinh  kdt 
and 

aH  sinh  k(dt+y) 

vw  *  T  (im sin(kx-at)  (4) 

where  a  =  wave  radian  frequency,  Hw  =  wave  height,  k  =  2 n/X,  \=  wavelength, 
dt  =  water  depth  and  x,y  are  horizontal,  vertical  coordinates  with  their  origin 
at  the  mean  water  level  directly  above  the  hinge. 

The  overturning  moment  acting  on  the  pile  as  a  result  of  (steady)  wind. 

My,  is  evaluated  using  a  drag  coefficient  approach.  An  approximate  expression 
for  wind  moment  can  therefore  be  written  in  the  form: 

\  ■  >/%2  «s2e  -  d2>  P,  c,  V,2  *  1/2  ib  COSS  pa  Cb  Ah  Ua2  (5) 

in  which  distances  d  and  are  shown  in  Fig.  2,  pa  =  air  density,  =  drag 
coefficient  of  the  pile  in  air,  s  drag  coefficient  of  daymark  boards,  Ab  = 
area  of  boards,  and  U  =  wind  velocity.  The  first  term  on  the  right  hand  side 

Q 

(RHS)  is  set  to  zero  should  the  pile  become  entirely  submerged,  and  the  second 
is  zero  when  the  boards  are  sacrificed. 

During  a  collision,  the  barge  contact  force  contributes  a  moment  about  the 
hinge,  Mg,  which  is  of  the  form 

"a  *  pb  *c  (6) 

where  Fg  =  barge  contact  force,  =  distance  from  hinge  to  point  of  contact  and 
6C  =  angle  between  pile  direction  and  direction  of  barge  force. 

The  general  pile  dynamic  expressions  given  by  Eqs.  1-6  serve  as  the  basis 
for  modeling  the  specific  conditions  stated  in  the  design  criteria.  Equation 
specialization,  solution  approaches  and  computer  programs  based  on  the  mathe¬ 
matical  theory  are  discussed  in  the  following  subsections  for  each  application. 


PRELIMINARY  ANALYSIS 


Two  computer  models  were  developed  to  assist  in  the  Initial  determination  of 
CTPS  design  parameters  and  to  provide  a  preliminary  assessment  of  the  system's 
static  and  dynamic  characteristics.  The  computer  program  PILESTIFF  calculates  the 
initial  hinge  moment  stiffness  k-j  (see  Fig.  3)  necessary  to  meet  a  specified 
verticality  requirement  under  static  equilibrium  conditions.  The  program  PILEFREQ 
computes  a  specified  pile  system's  undamped  natural  frequency. 

PILESTIFF  is  helpful  in  the  early  stages  of  design  when  a  trial  value  for 
kj  is  needed.  This  can  be  obtained  by  ignoring  the  oscillations  induced  by 
waves  and  solving  the  corresponding  static  equilibrium  problem.  Under  static 
conditions  the  left  hand  side  (LHS)  of  Eq.  1  is  zero,  while  the  RHS  includes 
moment  contributions  given  by  Fig.  3  and  Eqs.  2,3  and  5  (no  barge  contact).  The 
wave  fluid  velocity  contribution,  given  by  Eqs.  4  to  the  relative  velocity  (ur) 
in  Eq.  3  is,  however,  zero.  For  small  angles,  the  resulting  moment  equilibrium 
equation  is  easily  rearranged  to  provide  the  following  formula  for  k^: 


k!  *  V*  *  hwP  *  *  ko.Wi2  * 


(7) 


The  user  of  PILESTIFF  specifies  the  static  angle  desired,  all  other  CTPS  design 
parameters,  and  the  wind,  current  and  weight  loading.  PILESTIFF  then  uses  Eq.  7 
to  compute  the  value  for  k-j . 

PILEFREQ  calculates  the  CTPS  undamped  natural  frequency  which  is  useful  in 
identifying  potential  resonant  situations  from  wave  loading  or  other  sources  of 
periodic  forcing.  The  angular  momentum  equation  is  used  in  which  wave  excitation 
and  fluid  damping  (Eq.  3),  wind  forcing  (Eq.  5),  and  barge  contact  (Eq.  6)  are 
neglected  from  the  RHS.  The  remaining  hinge  moment  (Fig.  3)  and  weight  moment 
(Eq.  2)  terms  are  linearized  yielding  the  harmonic  oscillator  equation  from  which 
the  natural  frequency  u»Q  is  easily  identified  as 
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(8) 


wo  =  ^kr£mW2  ‘  zVV^HT-1 
where  IHT  =  gmp2p2  +  \Z  n>2  +  Pwcmdp2d3»  mp  =  mass  of  Pile  and  =  mass  of 
load. 

PILEFREQ  computes  natural  frequency  (and  period)  using  Eq.  8  and  CTPS 
design  parameters  supplied  by  the  user. 

OPERATING  CONDITIONS 

The  computer  program  OPPILE  was  developed  to  model  pile  dynamics  during 
operating  conditions.  The  principal  use  of  this  model  is  to  determine  whether 
the  CTPS  designs  under  consideration  meet  the  verticality  requirement  specified 
in  Table  2.  Thus  the  program  predicts  inclination  angle  and  hinge  moment  response 
for  specified  pile  dimensions,  hinge  stiffness  and  load  conditions. 

The  model  is  based  on  the  angular  momentum  equation,  Eq.  1,  with  piecewise 
linear  hinge  stiffness  behavior  (see  Fig.  3)  and  loads  due  to  weight,  current, 
wave,  pile  motion  relative  to  the  fluid,  and  wind  as  provided  by  Eqs.  2,3  and  5. 
Since  the  operating  restrictions  require  that  the  pile  be  nearly  vertical,  small 
angle  approximations  are  used  and  the  boards  are  assumed  not  to  submerge. 

Coupling  with  lateral  motion  is  neglected.  The  possibility  of  transverse 
excitation  by  vortex  shedding  was  investigated,  and  motion  due  to  this  source 
was  found  negligible  for  full  scale  CTPS's  meeting  the  design  criteria.  Reynolds 
numbers  (MixlO^)  exceed  the  critical  Reynolds  number  so  that  the  wake  is  fully 
turbulent  with  little  coherent  vortex  street  structure  as  discussed  by  Weigel 
(1964)  and  others.  Lift  coefficients  (C^)  in  this  range,  as  reported  for  example  by 
Sarpkaya  (1976)  and  Schewe  (1982),  drop  to  the  range  .03  <  <  .20.  In  addition 

the  natural  frequency  of  pile  oscillation  was  found  to  be  much  slower  than  the 
shedding  frequency,  so  the  system's  dynamic  response  is  very  small. 

The  specialized  dynamic  equations  are  solved  using  a  Runge-Kutta  numerical 
technique.  At  each  time  step,  the  along-pile  integrations  required  by  Eq.  3  are 
completed  using  the  trapezoidal  rule.  No  stability  problems  were  encountered, 


and  accurate  results  were  obtained  for  time  steps  less  than  a  tenth  of  the  wave 
period. 

The  program  user  must  supply  pile  and  hinge  design  parameters,  wind  and 
current  velocities,  water  depth,  and  wave  height,  period  and  length.  A  utility 
program,  WAVELENGTH,  was  written  to  assist  the  user  in  specifying  consistent  wave 
parameters.  Specifically  the  program  calculates  wavelength  for  user  specified 
depth  and  wave  period  by  solving  the  following  transcendental  equation  from 
small  amplitude  wave  theory: 

nT2 

A  =  Jj-  tanh  (2ndt/\)  (9) 

where  A  =  wavelength,  g  =  gravitational  acceleration  and  T  =  wave  period. 
Wavelength  A  is  computed  using  the  Newton-Raphson  iteration  method. 

When  input  to  OPPILE  is  complete,  the  program  calculates  and  prints  out 
time  series  for  inclination  angle  0  and  hinge  moment  M^.  Steady  state  response 
is  generally  achieved  within  3  wave  periods. 

HURRICANE  CONDITIONS 

The  computer  program  HURPILE  was  developed  to  model  pile  dynamics  during 
hurricane  conditions  such  as  those  described  in  Table  3  of  the  design  criteria. 
HURPILE  is  actually  very  similar  to  OPPILE.  The  main  difference  is  that  HURPILE 
is  not  limited  to  small  inclination  angles.  The  pile  may  be  entirely  submerged 
and  rotate  up  to  0  =  90  deg  without  loss  of  accuracy.  The  large  angle  capability, 
however,  necessitates  approximately  three  times  the  computer  time.  Another 
important  characteristic  is  that  the  boards  are  assumed  to  have  been  sacrificed. 
From  the  user's  point  of  view,  input  and  output  format  are  virtually  identical 
to  that  of  OPPILE. 
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COLLISION  CONDITIONS 

The  computer  programs  COLPILE  and  RECPILE  were  developed  to  model  pile 
dynamics  during  collision  conditions.  COLPILE  was  used  to  predict  angular 
position,  barge  loads  and  hinge  reactions  during  a  head-on  (worst  case)  barge 
collision.  RECPILE  predicts  pile  motion  and  hinge  moment  as  the  pile  recovers 
to  an  upright  position  after  the  CTPS  has  been  overrun.  A  discussion  of  COLPILE 
appears  immediately  below  and  is  subsequently  followed  by  a  description  of  RECPILE. 

The  complete  collision,  as  modeled  by  COLPILE,  consists  of  a  sequence  of 
processes.  Initially  there  is  impact  of  the  top  of  the  barge  bow  with  the  pile, 
then  sliding  of  the  pile  occurs  along  the  top  of  the  bow  and  the  bow  face.  Next 
there  is  impact  of  the  bottom  of  the  bow  rake  with  the  pile  followed  by  sliding 
along  this  point.  Lastly,  the  tip  of  the  pile  slides  along  the  barge  bottom  before 
being  released.  The  program  analyzes  these  processes  in  chronological  order. 

The  major  assumption  throughout  the  collision  analysis  is  that  because  the 
barge  is  so  massive,  it's  motion  is  essentially  unaffected  by  the  collision. 

Barge  speed  therefore  remains  constant.  In  addition,  it  is  assumed  that  the 
boards  are  sacrificed,  and  wind  and  wave  forces  are  considered  negligible  in 
comparison  with  collision  forces. 

The  constant  barge  speed  condition  enables  the  pile  kinematics  to  be  analyzed 
independently  of  the  forces  involved.  Since  the  horizontal  velocity  component  of 
the  pile  contact  point  must  equal  the  barge  speed,  0,  6  and  8  may  be  determined 
as  function  of  time  from  the  problem  geometry. 

Next,  CTPS  dynamics  are  analyzed  using  the  rate  of  change  of  angular  momentum 
equation,  Eq.  1,  in  which  the  moment  sum  includes  the  hinge  moment  M^,  the 
gravitational  moment  Mg,  the  fluid  force  moment  due  to  current  and  relative  pile 
motion  (no  waves)  Mg,  and  the  barge  moment  Mg.  Using  the  kinematic  results,  M^, 

Mg  and  Mg  are  evaluated  using  Fig.  3,  Eq.  2  and  Eq.  3,  respectively.  Similarly 
the  rate  of  change  of  angular  momentum  term,  1^8,  in  Eq.  1  is  calculated  from  the 
barge  kinematics.  Eq.  1  is  then  used  to  compute  barge  moment  MR. 


Note  that  unlike  OPPILE  and  HURPILE,  Eq.  1  is  not  “solved"  in  the  usual 
sense  of  evaluating  0  =  0(t).  The  constant  barge  speed  constraint  and  geometry 
determines  angular  position  independently,  and  Eq.  1  is  simply  used  to  calculate 
the  unknown  barge  moment,  Mg,  term. 

Using  the  Mg  result  and  colli  son  geometry,  Eq.  6  is  applied  to  compute  the 
barge  contact  force  Fg.  Having  determined  Fg,  the  linear  momentum  equations 
(Newton's  Second  Law)  are  applied  vertically  and  horizontally  to  evaluate  the 
hinge  reaction  forces  Ry  and  R^,  respectively. 

While  the  collision  analysis  described  above  is  theoretically  correct 
throughout  the  sequence  of  collision  processes,  it  is  convenient  to  modify  the 
approach  somewhat  at  the  two  instants  of  impact.  At  these  times,  the  weaker  non- 
impulsive  loads  due  to  fluid  motion  and  the  hinge  moment  are  neglected,  and 
COLPILE  uses  the  impulse-momentum  form  of  the  reduced  equations  of  motion.  Eq.  1, 
for  example,  becomes 

t. 

(IHe)t  -  (^t.  =  2/  1  M  dt  (10) 

f  i  t  • 

where  t.  and  t^  are  initial  and  final  times,  respectively,  bracketing  the  impact 
process  and  the  RHS  includes  only  impulsive  moments.  The  same  steps  previously 
outlined  are  taken  yielding  results  for  barge  moment  impulses  and  barge  and  hinge 
reaction  force  impulses. 

The  program  user  must  supply  COLPILE  with  CTPS  design  parameters,  barge 
dimensions  and  speed,  current  velocity  and  water  depth.  The  user  also  specifies 
the  interval  between  times  for  which  output  is  desired.  Results  consist  of  angle, 
hinge  moment,  barge  moment  (or  moment  impulse),  barge  force  (or  force  impulse), 
and  reaction  forces  (or  force  impulses). 

The  program  RECPILE  is  used  to  model  pile  recovery  after  being  released 
from  beneath  the  barge  bottom.  Initial  position  is  specified  by  the  user  (from 
COLPILE  results),  then  RECPILE  calculates  angular  motion  and  hinge  moment  as  the 

pile  returns  to  the  upright  position.  RECPILE  is  actually  a  modification  of 
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HURPILE  thus  making  use  of  HURPILE's  large  angle  capability.  Changes  include 
omitting  wind  and  wave  excitation  and  allowing  the  user  to  specify  the  initial 
conditions  for  0. 

The  user  of  RECPILE  must  input  the  CTPS  design  parameters,  current  and 
water  depth  as  well  as  the  initial  .nclination  angle.  The  program  responds  by 
calculating  and  printing  out  time  series  for  9  and  M^. 


III.  HINGE  COMPONENT  DEVELOPMENT 


INITIAL  CONCEPT  DEVELOPMENT 


At  first,  the  very  simple  and  inexpensive  rubber  tube  hinge  suggested  by 
Miller  (1982)  was  considered.  In  this  concept,  flexibility  is  achieved  by  using 
a  section  of  rubber  tube  as  shown  in  Fig.  4.  Piecewise  linear  behavior  is  obtained 
by  having  the  wall  thickness  thin  enough  so  that  the  compressive  side  buckles  at 
large  angles. 

The  moment  as  a  function  of  angle  relationship  for  the  rubber  section  at  small 
angles  can  be  determined  from  beam  theory.  This  approach  can  therefore  be  used  to 
estimate  the  initial  stiffness  k^.  For  small  angles, 


Uh  _  k  _  El 
6  “  *1  "  i 

r 


(H) 


=  (E/Ar)[n/4(rJ-r})], 


where  E  *  Young's  modulus,  I  =  area  moment  of  inertia  for  the  rubber  tube  cross- 
section,  and  £r,  rQ  and  r.  are  the  tube  dimensions  shown  on  Fig.  4. 

If  the  pile  consists  of  a  12  inch  wooden  pole,  the  computer  models  indicate 
that  a  k-j  value  of  at  least  3.5  x  10^  ft-lbs/rad  is  necessary  to  meet  the  operating 
inclination  limit  of  the  design  criteria.  Using  r.  =  6  in,  rQ  =  9  in,  £r  =  18  in 
(commensurate  with  a  relatively  thick-walled  section)  and  E  =  1200  psi  (which  is 

4 

an  extremely  stiff  example),  k^  is  estimated  to  be  2.3  x  10  ft-lbs/rad.  Since 
this  result  is  less  than  an  order  of  magnitude  smaller  than  the  necessary  value, 
the  rubber  hinge  concept  was  not  developed  further. 

It  became  clear  that  achieving  sufficient  initial  stiffness  would  be  a  major 
obstacle.  To  address  this  problem,  other  concepts  involving  a  very  large,  central 
spring  were  considered  as  discussed  by  Swift  and  Baldwin  (1984).  The  central 
spring  would  be  housed  within  the  pile  itself.  Thus  the  next  phase  of  modeling 
activity  focused  on  the  analysis  of  systems  using  18  inch  steel  pipe  for  the  pile 
section. 
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PRELIMINARY  MOOELING 


The  computer  models  were  applied  to  a  trial  case,  “Design  CTPS"  in  order  to 
establish  guidelines  for  further  hinge  component  selection  and  development. 

A  sumary  of  the  tests  and  conclusions  is  presented  here,  while  details  are  in¬ 
cluded  in  the  interim  report  by  Swift  and  Baldwin  (1984).  The  "Design  CTPS" 
consisted  of  an  18  inch.  Schedule  20  steel  pipe  pile  with  lengths  chosen  for  the 
maximum  operating  water  depth  of  30  ft.  Wind  velocity,  current  velocity  and  wave 
heights  were  set  at  maximum  values,  and  pile  performance  was  computed  for  all 
design  conditions  and  some  typical  situations  not  explicitly  addressed  by  the 
design  criteria. 

An  initial  stiffness,  ,  value  of  approximately  6x10^  ft-lbs/rad  (or, 
equivalently,  10^  ft-lbs  at  10  deg)  was  found  necessary  to  meet  the  maximum 
inclination  limit,  under  worst  case  operating  conditions,  of  10  deg.  This  k^ 
value  results  in  a  5  deg  angle  under  the  corresponding  static  environment,  that  is, 
the  same  maximum  current  and  wind  but  no  waves.  The  optimum  breakpoint  angle  was 
found  to  be  10  deg  as  this  provides  the  largest  stiffness  during  operating  con¬ 
ditions  and  immediate  stiffness  reduction  beyond  the  operating  range.  A  large 
angle  stiffness,  k£,  of  approximately  1/10  k^  resulted  in  prompt  pile  recovery  after 
collision  and  could  be  reduced  further.  Since  a  minimum  moment  of  approximately 
10^  ft-lbs  is  necessary  to  initiate  recovery  at  90  deg,  ^  should  not  normally 
be  negative. 

Impact  of  the  barge  with  the  pile  at  the  bottom  of  the  bow  rake  was  found 
to  be  very  sensitive  to  bottom  clearance.  When  hinge  clearance  was  less  than 
approximately  7  ft,  force  impulses  on  the  pile  and  reaction  forces  supporting 
the  hinge  base  became  excessive.  Thus  the  mudline  clearance  of  3  ft  specified  in 
Table  4  could  result  in  impact  forces  causing  local  damage  at  points  of  contact 
and  possible  loosening  of  the  base.  Typical  barges  using  channels  in  the  Houston 
area,  on  the  other  hand,  have  drafts  of  less  than  12  ft.  For  this  more  common 
situation,  collision  loads  were  determined  to  be  tolerable. 


FINAL  CONCEPT  DEVELOPMENT 

The  preliminary  computer  analysis  had  shown  that  the  hinge  component  must 
meet  certain  performance  requirements  in  order  for  the  entire  CTPS  to  meet  the 
overall  design  specifications.  For  the  CTPS  to  meet  the  minimum  verticality 
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requirement,  the  initial  stiffness  constant,  k-j ,  should  be  approximately  6x10 
ft-lbs/rad,  and  the  breakpoint  angle  is  best  at  10  deg.  The  large  angle  stiffness 
is  less  critical  and  may  be  roughly  between  zero  and  1/10  k-j.  Lastly,  the  hinge 
must  be  able  to  rotate  down  from  the  vertical  a  full  90  deg  for  all  horizontal 
angles,  and  the  attachment  must  be  able  to  sustain  large  impact  loads. 

Several  hinge  concepts  were  considered,  and  their  potential  performance 
characteristics  were  compared  with  these  hinge  criteria.  (Examples  of  hinge 
component  designs  considered  but  not  adopted,  besides  the  rubber  tube  hinge,  are 
reviewed  in  Appendix  8).  The  system  found  to  be  best  able  to  meet  the  hinge 
criteria  is  the  central  universal  joint,  prestressed  peripheral  stay  concept 
illustrated  in  the  Fig.  5  schematic.  The  two  axes  of  the  universal  joint  are 
off-set  vertically  to  prevent  binding  at  any  point  in  the  range  of  articulation. 
Four  stays  are  attached  to  the  base,  led  upwards  over  spreaders  and  into  the 
(hollow)  pile,  and  connected  to  a  central  spring.  The  spring  is  prestressed  so 
that  all  stays  are  equally  taut  when  the  pile  is  vertical.  The  arrangement  is 
such  that  if  there  is  an  angular  change  from  the  vertical,  the  stay  on  the  outside 
of  the  bend  immediately  takes  up  the  entire  prestress  force.  Thus  a  large  re¬ 
storing  moment  is  generated  at  very  small  angles. 

The  prestress  moment  is  essentially  rgFs  where  rs  =  radius  of  stay  location 
and  F$  is  the  prestress  force.  Ideally  this  is  imposed  for  the  smallest  deviation 
from  the  vertical.  Because  the  real  system  is  non-ideal  (due  to  friction  and 
compliance  of  the  components),  some  finite  angle  change  must  occur  before  MH 
increases  to  rgFg.  To  be  sure  that  the  necessary  initial  stiffness  is  achieved, 
the  prestress  moment  is  set  at  the  computer  model  hinge  moment  evaluated  at  the 


Central 
Spring  Force 


limit  of  the  operating  range  (that  is,  10  deg).  The  prestress  force  is  therefore 
designed  to  be 

Fs  =  Vrs  (Hh(6  =  10°)) 

=  l/rs  k7  (lOn/180)  (12) 

Half-circular  sheaves  are  attached  below  in  order  to  maintain  stay  moment 
arm  with  respect  to  the  hinge  axis.  With  the  arm  distance  fixed,  the  hinge  moment 
normally  increases  slowly  at  large  angles  as  the  prestressed  spring  is  extended 
further. 

The  "spring"  used  must  have  the  capacity  for  providing  the  prestress  force 
given  by  Eq.  12  and  have  the  elastic  range  necessary  for  hinge  angles  up  to  90  deg. 
The  "spring"  must  also  fit  inside  the  pile.  Possible  options  include  rubber  bands, 
an  array  of  nylon  ropes  or  a  piston/cylinder  (air  bag  sealed)  concept.  Further 
discussion  of  "spring"  options  is  included  in  Appendix  C. 


IV.  CTPS  DESIGN  AND  MODEL  CONSTRUCTION 


PILE  SYSTEM  DESIGN 

The  central  universal  joint/peripheral  stay  hinge  concept  was  incorporated 
into  the  pile  system  design  shown  in  Fig.  6.  CTPS  parameters  were  selected  so 
that  the  system  would  be  operable  in  30  ft.  of  water  as  specified  in  the  design 
criteria.  The  pile  supports  two  angled  6  ft  by  6  ft  daymark  boards  (face  view 
of  one  shown  in  Fig.  6)  and  a  load  consisting  of  a  light/battery  pack.  These 
are  located  7  ft  above  the  high  water  level  as  recommended  in  the  Coast  Guard  aid- 
to-navigation  manual. 

The  pile  itself  consists  of  an  18  inch  steel  pipe  and  contains  the  spring 
(one  of  the  options  discussed  previously).  The  spring  is  attached  at  the  lower 
end  to  the  stays  just  above  a  central  stay  guide.  The  upper  end  is  secured 
to  a  worm  gear  arrangement  for  prestressing  the  system  and  adjusting  individual 
stay  tautness.  Thus  the  stay  controls  are  accessible  from  the  top. 

The  hinge  design  follows  the  central  universal  joint/prestressed  peripheral 
stay  concept  described  in  the  previous  section.  The  spring  pretension  force  was 
evaluated  using  computer  model  results  and  Eq.  12.  Our  preliminary  computer 
analysis  indicated  that  if  k^  is  determined  using  a  5  deg  inclination  angle 
criteria  under  static  worst  case  operating  conditions,  the  dynamic  response  should 
not  exceed  the  10  deg  maximum  angle  when  waves  are  present.  PILESTIFF  was, 
therefore,  used  to  compute  an  initial  stiffness  k^  =  577,770  ft-lbs/rad  from 
which  a  spring  prestress  force  F$  =  67,200  lbs  was  calculated  using  Eq.  12.  Be¬ 
cause  the  sheaves  maintain  stay  moment  arm  at  large  inclination  angles,  k2, 
though  small,  should  not  be  negative.  The  hinge  is  secured  to  a  heavy  base 
connected  to  an  embedded  pile  section. 

Since  this  design  was  developed  in  accordance  with  the  design  guidelines 
established  in  the  previous  section,  it  was  anticipated  that  the  system  would 
meet  the  design  criteria.  Further  evidence  was  obtained  by  testing  a  physical 

scale  model  and  conducting  a  complete  series  of  computer  simulations  as  described 
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1 .5  ft. 

Pig.  6.  The  CTPS  design  for  a  maximum  operating  depth  of  30  ft.  System 

weight  for  all  components  above  the  lower  shaft  is  4634  lbs.  Weight 
for  the  18  In  pile  only  Is  2890  lbs  ,  and  "load"  weight 
(including  all  hardware  as  well  as  the  beacon)  is  1744  lbs  and  is 
located  8.84  ft  above  the  lower  shaft.  For  all  components  above 
the  upper  shaft,  total  weight  is  3740  lbs  ,  pile  weight  is  2716  lbs 
"load"  weight  Is  1024  lbs  and  is  located  12.07  ft  above  the 
upper  shaft.  ?7 


in  following  sections.  Scaling  considerations  and  construction  details  for  the 
physical  model  are  presented  in  the  remainder  of  this  section. 

SCALING  CONSIDERATIONS 

Planning  for  the  construction  and  testing  of  the  physical  model  required  that 
scaling  decisions  be  made.  The  model  would,  of  course,  be  built  so  that  its 
shape  is  scaled  geometrically  and  mass  and  weight  are  proportional  to  Lp^  where 

Lf  =  (^ModeV^p^ull  Sized'  A  more  comP1ete  scaling  system,  however,  needed 
to  be  devised  in  order  to  establish  other  parameters  such  as  speeds  and  stay 
tensions. 

In  view  of  the  fact  that  both  inertial  and  gravitational  forces  play  a 
crucial  role  in  the  pile  dynamics,  it  was  decided  that  Froude  scaling  would  be 
adopted  for  the  water  channel  and  pool  experiments.  Holding  Froude  number 

(»  U/Jgd^  =  (inertial  forces/gravitational  forces)^2)  constant  in  a  constant 
gravitational  environment  results  in  velocities  and  time  scales  proportional  to 

JTr.  Consequently  IH  8  or(length)2(mass)(time) "2  cr  LpA;  from  Eq.  2,  WQ  or  Lp4, 
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and  from  Eq.  3,  Mc  a  Lp  ,  . 

In  order  that  the  dynamic  scaling  of  Eq.  1  be  consistent,  it  is  also 
4  3 

required  that  Mu  or  L„  and  therefore  F.  or  L_  .  It  should  be  noted  that  this 
Hr  s  r 

represents  a  reduction  at  the  model  size  from  that  obtained  by  straight  geometric 

scaling.  Since  In  general,  Fg  a  (elastic  constant)(strain)(cross-section  area), 

using  the  same  material  and  geometrically  Imposing  the  same  initial  strain 
2 

causes  F  or  L  .  Thus  if  spring  material  is  unchanged,  F  must  be  further  reduced 
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to  comply  with  F$  a  Lp  scaling.  This  may  be  achieved  by  reductions  in  initial 
strain,  cross-section  area  or  both. 

The  system  of  scaling  outlined  here  was  used  as  the  basis  for  specifying 
parameters  for  the  physical  scale  model  fabrication  described  in  the  next  subsection 


The  scaling  method  was  also  employed  in  all  physical  model  experiments  as 
discussed  in  the  following  sections. 

MOOEL  CONSTRUCTION 

A  physical  scale  model  was  constructed  to  correspond  to  the  full  scale  design 
shown  in  Fig.  6.  A  scale  ratio  Lr  =  1/15  was  chosen  so  that  the  specified  high 
water  operating  condition  (depth  =  30  ft)  could  be  tested  at  model  scale  in  the 
Coast  Guard  water  channel  (depth  =  2  ft).  Construction  drawings  of  the  model 
are  provided  in  Figs.  7-10,  and  a  photograph  of  the  hinge  assembly  is  shown  in 
Fig.  11. 

In  order  to  minimize  complicated  fabrication  tasks,  make  use  of  standard 
sections  and  to  provide  sufficient  local  strength,  major  components  of  the  model 
are  overly  thick.  Consequently  material  substitutions  were  necessary  to 
maintain  correctly  scaled  weights.  The  pile  itself  is  made  of  PVC  pipe  rather 
than  steel,  and  aluminum  is  substituted  for  steel  in  the  hinge  assembly. 

A  very  simple,  conveniently  available  rubber  band  is  used  as  the  "spring." 

This  is  pre-stressed  to  a  force  of  19.9  lbs  (=  lr3(Fs)Fuli  Scale  =  67,200) 

corresponding  to  a  k1  value  of  11.4  ft-lb/rad  (=  Lr^(k1)pui1  Scale  ®  0/15)4  577,700). 
A  reel  for  prestressing  the  rubber  band  is  positioned  at  the  top.  Since  stay 
moment  arm  with  respect  to  the  hinge  axis  is  maintained  by  the  sheaves  at  large 
angles  and  spring  tension  increases  slightly  with  angle,  k£  is  expected  to  be  small 
but  positive.  Its  exact  value,  however,  will  depend  greatly  on  internal  friction 
and  is  difficult  to  quantify  theoretically. 

Fine  adjustment  of  stay  tension  is  provided  in  the  model  by  use  of  turnbuckles 
located  at  the  base.  Each  stay  is  led  through  a  base  attachment  eye  then  to  a 
turnbuckle  which  lies  horizontally  just  above  the  base.  This  system  enables 
convenient  tuning  for  model  experiments.  The  full-scale  design,  of  course,  will 
have  all  tension  gear  accessible  from  the  top  as  noted  previously. 

Flexible  wire  rope  was  chosen  for  the  stays  because  of  its  strength  and 
availability  at  an  appropriately  small  diameter.  Kevlar  is  recommended  for  the 


Hinge  component  parts  used  In  the  upper  shaft  assembly 
Drawings  are  full  scale,  and  dimensions  are  in  inches. 
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Fig.  11.  Physical  model  hinge  assembly.  A  vertical 
scale  is  provided  which  has  units  of  inches 
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V.  HINGE  MOMENT/ANGLE  EXPERIMENTS 


OBJECTIVE 

With  the  properly  scaled  pile  model  built  the  next  step  was  to  obtain  data 
to  determine  the  actual  hinge  moment  (M^)  vs.  angle  (0)  relationship.  The 
testing  produced  plots  of  M^  vs.  9  which  defines  the  hinge  behavior.  Values  of 
k-j  and  kg  were  evaluated  and  compared  with  those  obtained  from  the  theoretical 
and  design  work. 

INSTRUMENTATION  AND  PROCEDURES 

Prior  to  the  actual  measurement  of  the  M^  vs.  0  relationship  the  spring  was 
prestressed  to  the  design  value  of  Fs  =  19.9  lbs.  The  actual  main  spring  element 
incorporated  in  the  model  was  pulled  in  an  Instron  testing  machine  to  generate 
an  axial  force/elongation  plot.  To  obtain  the  desired  preload  in  the  model,  the 
spring  element  was  elongated  to  the  value  corresponding  to  Fg  =  19.9  lbs. 

With  the  spring  preloaded,  a  sequence  of  angular  displacement,  perpendicular 
force  and  moment  arm  measurements  were  recorded.  This  data  was  subsequently 
reduced  to  provide  MH  vs.  0  plots  and  values  of  and  kg. 

The  actual  procedure  and  equipment  employed  in  the  measurements  are  shown 
schematically  in  Fig.  12.  The  pile  was  secured  to  a  rigid  table.  An  angle 
chart  was  oriented  such  that  the  vertex  of  the  angle  chart  was  in  line  with  the 
axis  about  which  the  measurement  was  to  occur.  A  pointer  was  attached  to  the 
pile  to  serve  as  an  angle  indicator.  The  force  was  applied  at  some  distance 
above  the  axis  of  rotation  by  means  of  a  block  and  line  arrangement  as  shown  in 
Fig.  12.  The  load  cell  in  the  line  was  used  to  measure  the  force.  Adjustments 
were  made  in  the  line-block  system  to  insure  that  the  force  was  being  applied 
perpendicular  to  the  pile  at  each  angle  0  where  data  was  recorded.  The  load 
cell  was  an  Interface  super  mini  load  cell  with  a  range  from  0  to  10  lbs  tension 


or  compression.  It  was  powered  with  a  5  volt  dc  power  supply  and  the  output 
was  indicated  on  a  Fluke  digital  multimeter.  Data  was  obtained  at  2°  increments 
until  10°  and  then  at  15°  increments  until  the  geometric  limits  of  the  measurement 
system  were  achieved,  typically  60°. 

Three  pieces  of  information  were  obtained  to  provide  one  point  on  the  vs. 

0  plots  -  angle,  moment  arm  and  force.  The  moment  arm  times  the  force  provided  the 
measurement  Mm.  This  value  had  to  be  corrected  for  the  gravitational  moment  Mq 
given  in  Eq.  2.  The  gravitational  moment  was  most  significant  at  large  values  of 
8.  The  correction  was  made  using  the  following  equation 

MH  =  Mm  +  MG'  <13> 

RESULTS 

This  procedure  was  employed  for  three  cases:  1)  rotation  about  the  lower 
axis,  2)  rotation  about  the  upper  axis  and  3)  an  oblique  case.  In  the  oblique 
case  the  pile  is  displaced  about  both  axes  simultaneously.  The  data  obtained 
for  each  case  was  reduced  and  plots  of  MH  vs.  0  for  each  of  the  cases  are  shown 
in  Figs.  13,14  and  15.  The  actual  data  points  are  indicated  on  each  plot  along 
with  a  least  squares  fit  of  a  piecewise  linear  curve  having  a  break  at  0  =  10  deg. 
Values  of  k-j  and  k2  were  calculated  from  this  least  squares  analysis  and  are 
provided  on  the  plots. 

The  values  of  k^  and  k2,  summarized  in  Table  5,  vary  from  case  to  case. 

The  point  to  note  is  that  the  values  of  k^  for  the  lower  axis  and  oblique  cases 
are  within  10%  of  the  design  values.  The  values  of  k^  for  the  upper  axis  case  is 
higher  and  therefore  conservative  and  acceptable.  Though  highly  variable,  all  k2 
values  correspond  to  hinge  moments  sufficient  to  return  the  pile  to  the  vertical 
position  from  a  horizontal  starting  point. 

The  difference  between  the  upper  and  lower  shafts  was  considered  an  artifact 
of  the  pretensioning  fine  tuning.  The  stays  are  adjusted  in  pairs  to  obtain 
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Fig.  13.  Hinge  moment  about  the  upper  axis  vs.  angle.  A  least 
squares  fit  of  the  data  yields  k,  »  16.06  ft-lb/rad  and 
k2  *  -.4352  ft-lb/rad  (or  k,  «  0.2803  ft-lb/deg  and 
-  kf  *  -0.0076  ft-lb/deg).  Tne  resulting  piecewise 
linear  curve  Is  shown. 
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Fig.  14.  Hinge  moment  about  the  lower  axis  vs.  angle.  A  least 
squares  fit  of  the  data  yields  k,  »  10.37  ft-lb/rad  and 
k-  -  0.0152  ft-lb/rad  (or  k,  ■  0.181  ft-lb/deg  and 
k-  *  0.0142  ft-lb/deg).  The  resulting  piecewise  linear 
carve  Is  shown. 
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Fig.  15.  Hinge  moment  oblique  to  the  hinge  axes  vs.  angle.  A  least 
squares  fit  of  the  data  yields  k.  *  10.46  ft-lb/rad  and 
■  .5141  ft-lb/rad  (or  k,  ■  0.1825  ft-lb/deg  and 
kg  ■  0.0084  ft-lb/deg).  Tne  resulting  piecewise  linear 
carve  Is  shown. 
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VI.  WATER  CHANNEL  EXPERIMENTS 


OBJECTIVES 

The  main  objective  was  to  determine  if  the  model  could  satisfy  the 
verticality  requirement  in  a  steady,  uniform  current.  Secondary,  yet  important 
information  regarding  the  hydrodynamic  moment  and  drag  coefficient  was  also  an 
objective  of  these  experiments.  The  uniform  current  testing  was  performed  at  the 
U.S.  Coast  Guard  Academy  Circulating  Water  Tunnel  Facility  (CWT)  in  New  London,  CT. 
Two  types  of  experiments  were  performed  to  acquire  the  desired  data.  In  the  first 
experiment  the  pile  inclination  angle  was  measured  as  a  function  of  the  current 
speed.  This  allowed  the  observation  and  quantification  of  the  verticality  con¬ 
dition.  In  the  second  experiment  the  hydrodynamic  moment  was  measured  as  a 
function  of  current  speed.  This  relationship  enabled  the  calculation  of  the  drag 
coefficient. 

INSTRUMENTATION  AND  PROCEDURES 

In  both  experiments  the  pile  base  was  secured  to  the  bottom  of  the  test 
section  of  the  CWT.  A  bridge  structure  which  straddled  the  test  section  was  put 
in  place  and  secured.  This  structure  provided  a  fixed  reference  frame  for  measur¬ 
ing  pile  angle  and  attaching  the  load  cell  for  moment  measurements  as  shown  in 
Figs.  16  and  17.  The  water  speed  (Uc)  was  varied  from  0.6  ft/sec  (2.3  ft/s  full 
scale)  to  1.4  ft/sec  (5.4  ft/s  full  scale).  The  measured  speed  was  determined  by 
timing  a  wooden  float  in  the  flow  over  a  10  foot  length  a  number  of  times  and 
calculating  an  average  speed. 

Angle  as  a  Function  of  Speed 

The  first  experiment  performed  was  the  angle  vs.  water  speed  test.  The  spring 
was  prestressed  and  fine  tuned  with  the  turnbuckles  to  achieve  pile  vertically, 
and  the  pile  base  was  secured  to  the  bottom  of  the  test  section.  The  horizontal 
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Fig.  17.  Circulating  Water  Channel  setup  for  hydrodynamic  moment 
vs.  speed  experiment.  Upper  schematic  shows  details 
of  load  cell  -  pile  arrangement. 
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distance  from  the  longitudinal  I-beam  bracket  (see  Fig.  16)  on  the  bridge 
structure  to  the  top  of  the  pile  was  recorded.  This  value  was  considered  a  re¬ 
ference  value  (dfef)  indicating  the  vertical  position  of  the  pile. 

The  water  speed  was  increased  to  a  value  in  the  range  listed  above  and  a 
steady  state  was  allowed  to  develop.  The  speed  was  then  measured  along  with  the 
horizontal  distance  (d^)  from  the  pile  top  to  the  bracket.  The  angle  was  de¬ 
termined  by: 


=  tan  9 
£P 


(14) 


where  d^  and  dre^  are  as  defined  above  and  Z is  the  distance  from  the  axis  of 
rotation  to  the  top  of  the  pile  (as  shown  in  Fig.  16).  This  procedure  was 
repeated  at  different  speeds  until  the  upper  limit  of  the  water  speed  range  was 
achieved. 


Hydrodynamic  Moment  Measurements 

In  the  second  experiment  the  hydrodynamic  moment  (Mc)  due  to  was  measured. 
The  pile  base  plate  remained  secured  to  the  CWT  test  section  bottom.  The  fore 
and  aft  stays  were  disconnected  allowing  unrestrained  pile  rotation  in  the  direction 
of  the  flow.  A  load  cell  was  bolted  to  the  bridge  bracket,  as  shown  in  Fig.  17, 
fixing  one  end  of  the  load  cell.  The  action  end  of  the  load  cell  was  then  connected 
to  the  top  of  the  pile  with  a  piece  of  light  nylon  line.  The  pile  was  moved  to  a 
vertical  position  and  the  nylon  line  secured.  The  location  of  the  line  on  the 
pile  relative  to  the  axis  of  rotation  was  recorded.  The  water  speed  was  then  in¬ 
creased  as  in  the  previous  experiment  and  the  value  of  the  force  was  recorded  for 
each  Ug.  The  load  cell  used  in  this  experiment  was  the  same  one  used  in  the  bench 
testing.  The  distance  Z  times  the  force  yielded  the  hydrodynamic  moment  (Mr) 
due  to  the  uniform  current  U-. 


RESULTS 


The  results  from  the  angle  vs.  current  experiment  are  plotted  in  Fig.  18 
where  current  speed  is  provided  for  both  the  model  scale  and  full  scale.  The 
main  point  to  note  here  is  the  low  value  of  the  angle  at  maximum  Ug.  The 
angular  displacement  of  the  pile  is  less  than  1  deg  from  the  vertical  at  the 
maximum  current.  This  is  well  within  the  verticality  limit  specified  in  the 
design  criteria.  From  this  perspective  the  experiment  was  highly  successful  and 
demonstrated  the  integrity  of  the  design  concept  of  high  initial  stiffness  to 
maintain  verticality  under  environmental  loading. 

The  computer  model  PILESTIFF  was  used  to  estimate  an  effective  stiffness  for 
the  test  conditions.  It  was  found  that  the  effective  during  the  trials 
averaged  about  twice  the  design  value.  Thus  in-water  performance  was  found  to 
be  better  than  both  theoretical  predictions  and  bench  tests. 

Hydrodynamic  moment  measurements  were  used  to  estimate  the  drag  coefficient, 
Cw,  and  to  compare  this  with  the  Cw  values  assumed  in  the  initial  computer 
modeling.  The  drag  coefficient  was  determined  using  the  following  equation. 


r  s  _ _ n. 

w  1.  ,,2, 


(15) 


Wp 

Mc(d/2) 


■1 


kuc<dpd) 

and  results  are  plotted  in  Fig.  19.  There  is  some  variability  in  the  Cw  values 
determined,  but  all  points  (with  average  =  0.886)  are  less  than  the  value  of 
Cw  =  1  used  in  the  computer  models. 

During  both  experiments  observations  were  made  regarding  the  lateral  motion 
of  the  pile.  Lateral  motion  is  considered  the  motion  perpendicular  to  the  flow 
direction  which  is  induced  by  vortex  shedding.  The  observed  lateral  motion  was 
very  limited,  approximately  1/16  inch  peak  to  peak  at  (Ur)__v  as  determined  by 
a  hand  held  scale. 


0.6  0.7  0.8  0.9  I .0  T.I  T.2  1.3  1.4 

^UC^Model  ^t/s^ 


^C^Full  Scale  (ft/s) 


'C'Full  Scale  VM,ul*s' 


Results  of  the  angle  vs.  speed  experiment. 
Maximum  angle  (8  max.)  corrected  for  wind 
using  OPPILE  (60  knots  scaled)  is  2.23  deg 
which  is  less  than  the  most  restrictive 
design  criteria. 


Ft g .  19.  Results  of  the  hydrodynamic  moment  vs.  speed  experiment 
Shown  are  the  hydrodynamic  moment  M-(+)  due  to  current 
and  the  drag  coefficient  C,,(o). 


The  results  of  these  experiments  were  encouraging.  The  ability  of  the  model 
to  provide  the  sufficient  initial  stiffness  under  the  maximum  uniform  current 
condition  indicated  the  concept  is  an  achievable  solution  to  the  design  problem. 
The  drag  coefficient  determined  was  conservative  with  respect  to  the  assumed 
value  used  in  the  theoretical  work.  Finally,  the  behavior  of  the  model  with 
respect  to  vortex  excitation  was  minimal. 


VII.  BARGE-PILE  COLLISION  EXPERIMENTS 


OBJECTIVES 

The  objectives  of  the  experiments,  performed  at  the  UNH  indoor  pool,  were 
twofold.  The  primary  objective  was  to  qualitatively  observe  and  record  the  impact 
processes  with  movies  and  still  photography.  The  important  point  to  observe  was 
the  return  of  the  pile  to  a  vertical  position  after  the  collision.  The  second 
objective  was  to  measure  the  impact  forces  thereby  quantifying  the  collision. 

This  required  instrumenting  the  barge  with  load  cells  and  recording  the  data  in  a 
permanent  form. 

FACILITIES,  INSTRUMENTATION  AND  PROCEDURES 

Facilities 

The  UNH  pool  is  75  ft.  long  with  an  approximate  depth  of  4  ft.  along  its 
shallow  edge.  Thus  a  platform  had  to  be  constructed.  When  the  pile  base  was 
attached  to  the  platform,  the  proper  amount  of  the  pile  (2  ft.)  was  below  the 
ambient  waterline.  The  platform  was  placed  approximately  50  feet  from  one 
end,  as  shown  in  Fig.  20,  to  allow  room  for  bringing  a  scale  model  barge  up  to 
speed  before  the  collision.  The  pile  was  secured  to  the  platform  with  the  upper 
hinge  axis  oriented  to  be  the  axis  of  rotation  during  the  collision  with  the  barge. 

The  barge  was  a  wooden  vessel  constructed  in-house.  Water  was  used  for 
ballast,  and  the  barge  was  partitioned  to  reduce  the  effect  of  sloshing.  When  the 
barge  was  entirely  filled  with  water,  the  (maximum)  draft  was  0.75  ft.  (11.25  ft 
full  scale).  The  forward  compartment  was  left  dry  for  housing  the  impact  in¬ 
strumentation. 

The  barge  was  pulled  through  the  water  by  means  of  a  hand  reel  to  wind  up 
a  tow  rope.  The  50  ft.  of  unobstructed  travel  was  sufficient  for  the  barge  to 
come  up  to  a  steady  speed  which  was  measured  by  timing  the  barge  as  it  passed  over 
the  platform. 
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Fig.  20.  Barge-pile  collision  experiment  set-up  at  the  UNH  pool 


Movies  and  still  photos  of  the  collision  process  were  made  both  above 
and  under  water.  A  movie  camera  enclosed  in  an  underwater  housing  was  located 
on  the  platform  directly  across  from  the  pile  hinge  with  the  film  plane  parallel 
to  the  plane  containing  the  barge  velocity  vector  and  the  pile.  The  camera  was 
also  mounted  above  water  at  an  angle  to  the  barge  direction  to  film  the  upper 
collision.  The  movie  camera  locations  are  designated  by  M  on  Fig.  20.  The 
photography  was  completed  with  a  set  of  still  exposures  taken  from  several 
locations  both  above  and  below  the  water. 


Impact  Force  Instrumentation 


Impact  instrumentation  included  force  sensors  consisting  of  two  cantilever 
beam  load  cells  which  were  fabricated  by  the  project  team.  Strain  gages  were 
attached  to  the  upper  and  lower  surfaces  to  measure  the  axial  strain.  The  load 
cells  were  essentially  two  active  arm  transducers  as  shown  in  Fig.  21.  The  bridge 
excitation  and  output  were  monitored  with  a  Validyne  CD- 19  carrier  demodulator. 

The  load  cells  were  connected  to  the  CD-19,  and  the  system  was  fine  tuned  and  cali 
brated. 


The  load  cells  were  mounted  inside  the  dry  compartment  in  the  barge  bow  area 
at  two  locations  (see  Fig.  22)  of  importance  in  the  collision  process.  One  was 
positioned  to  sense  the  upper  impact  (location  A),  and  the  other  was  set  to  sense 
the  lower  impact  (location  B).  Carriage  bolts  were  threaded  into  the  free  ends 
of  the  beams  and  allowed  to  extend  through  holes  to  beyond  the  outer  surface  of 
the  barge.  The  bolt  heads  were  covered  with  a  flexible  membrane  to  prevent 
flooding  of  the  compartment.  The  load  cells  were  oriented  in  a  manner  to  sense 
the  perpendicular  component  of  the  impact  force  as  indicated  in  Fig.  22. 

A  lightweight  aluminum  bumper  system  was  added  to  the  collision  zone  of  the 
barge  to  insure  that  the  Impact  was  recorded  as  It  would  be  difficult  to  hit  the 
bolt  head  alone  during  the  collision  experiments.  The  lightweight  bumper  system. 
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Fig.  21.  Schematic  of  Impact  force  measurement  system.  A  two 
active  arm  strain  gage  load  cell  Is  coupled  to  a 
signal  conditioning  device  (CO-19)  whose  output  Is 
displayed  and  stored  on  the  digital  storage  oscilloscope. 


Fig.  22.  Detail  view  of  the  barge  load  cell  -  bumper  system  for 
Impact  force  measurement.  Shows  the  pile-load  cell 
bumper  system  action  Is  perpendicular.  Location  A  Is 
the  upper  Impact  and  location  B  Is  the  lower  Impact. 


shown  in  Fig.  22,  consisted  of  a  pivot  point,  located  as  far  as  possible  away 
from  the  bolt  head  actuator,  an  aluminum  bumper  and  a  bracket  to  hold  the  free 
end  of  the  bumper  in  position  and  still  allow  movement  in  the  desired  direction. 

In  practice  the  angle  of  bumper  rotation  is  small,  and  the  point  of  impact  is 
very  close  to  the  bolt  head.  Thus  the  impact  force  on  the  bumper  was  essentially 
the  same  as  that  affecting  the  actuator- load  cell. 

Recording  the  impact  force  time  series  required  more  instrumentation.  A 
means  of  acquiring  the  signals  from  the  load  cells  had  to  be  configured.  As 
mentioned  earlier,  the  load  cells  were  calibrated  in  a  system  incorporating  a 
Validyne  CD-19.  The  recording  system  essentially  added  a  Nicolet  2096  digital 
oscilloscope  with  disk  storage  to  complete  the  measurement  system  as  shown  in 
Fig.  21.  The  electronics  was  set  on  a  cart  next  to  the  pool  side  and  connected 
to  the  load  cells  with  a  short  shielded  cable.  The  cart  was  moved  along  at  the 
same  rate  as  the  barge  allowing  a  minimum  amount  of  shielded  cable  to  be  used. 

The  oscilloscope  was  set  for  manual  trigger  and  was  triggered  just  before  the 
collision  thus  digitally  recording  the  impact  forces.  The  oscilloscope  memory 
was  subsequently  transformed  to  the  floppy  disk  creating  a  permanent  record  of 
the  event.  This  data  and  the  barge  velocity  allowed  for  a  quantification  analysis 
of  the  collision  process. 

RESULTS 

Observational 

Important  results  are  obtained  directly  from  observations.  The  experiments 
showed  that  there  is  no  physical  damage  to  the  model  during  impact  and  that  it 
returns  to  a  vertical  position  after  the  collision.  This  process  is  documented  on 
two  films  submitted  to  the  Coast  Guard  R&D  Center  at  Avery  Point,  CT.  The  still 
photographs  presented  here  in  Figs.  23  and  24  illustrate  the  four  steps  in  the 
collision  process:  1)  upper  impact,  2)  movement  of  the  pile  along  the  bow, 

3)  lower  impact  and  4)  the  movement  of  the  pile  along  the  barge  bottom  and  ultimate 


Plots  of  typical  upper  and  lower  impact  traces  are  shown  in  Pig.  25.  The 
upper  impact  time  series  is  truncated  at  approximately  7  lbs.  This  is  due  to 
the  bumper  coming  up  against  a  mounting  bolt  head  which  restricted  its  motion. 

A  sine  wave  was  fitted  to  the  remaining  correct  data  to  allow  the  peak  value  to 
be  estimated.  Knowing  the  peak  value  and  when  it  occurs  is  important  due  to 
the  nature  of  the  impact  process. 

The  impact  is  considered  to  have  started  at  t  =  0  and  continues  until  the 
time  when  the  force  reaches  its  peak  value.  This  represents  the  barge  hitting 
the  pile  over  some  finite  time  interval  during  which  the  load  cell  deflects  to 
its  maximum  excursion.  At  this  time,  the  pile  contact  point  has  accelerated  to 
the  barge  speed,  and  the  impact  is  over.  After  this  time  the  pile  commences  to 
move  away  from  the  barge,  and  the  load  cell  is  "released"  and  rings  until  the 
no  load  state  is  again  reached.  Thus  only  the  time  until  the  maximum  de¬ 
flection  is  considered  the  impact. 

From  the  time  series  of  the  impact  events,  the  integral  of  force  over  time 
(JFg(t)dt  =  force  impulse)  was  calculated  for  both  the  upper  and  lower  impacts. 
The  upper  impulse  is  0.124  lb-s  and  the  lower  impulse  is  0.306  lb-s.  COLPILE 
was  run  for  model  conditions  and  the  following  impulse  results  were  obtained: 
upper  impulse,  0.131  lb-s  and  lower  impulse,  0.236  lb-s.  These  results  are 
within  23%  at  the  largest  discrepancy,  indicating  that  the  COLPILE  program  is 
as  useful  as  a  design  tool. 


VIII.  COMPUTER  MOOEL  APPLICATION  TO  FULL-SCALE  PROTOTYPE 


PARAMETER  SELECTION 

The  computer  models  were  applied  to  the  CTPS  design  to  estimate  its  full 
scale  performance.  Linear  dimensions  for  computer  program  input  were  obtained 
directly  from  Fig.  6.  Other  model  parameters  were  specified  taking  into 
account  the  results  of  the  physical  scale  model  experiments.  Environmental 
conditions  simulated  -  operational,  hurricane  and  collision  -  were  selected  to 
conform  to  the  design  conditions  referred  to  in  Tables  2-4.  Input  to  the  models 
is  summarized  in  Table  6. 

The  hinge  initial  stiffness  coefficient  k-j  was  specified  as  the  design 
value  of  577,000  ft- lbs/rad,  and  the  breakpoint  angle  was  set  at  10  deg.  The 
bench  and  water  channel  tests  indicated  that  this  choice  of  stiffness  is  clearly 
achievable.  If  the  actual  k^  is  larger  than  the  computer  model  kp  which  the 
experiments  suggest  is  possible,  verticality  performance  will  be  improved.  Hence 
the  k-j  specified  is  a  "worst  case"  value.  When  the  actual  k-j  measured  in  the 
moment/angle  tests  approximated  the  computer  model  value,  k2  was  about  1/20  kp 
This  ratio  was  therefore  employed  in  specifying  k2  for  the  computer  simulations. 
Though  negative  k2's  were  observed  when  the  initial  stiffness  was  much  greater 
than  the  computer  model  value,  the  resulting  large  angle  hinge  moments  still  were 
approximately  the  same  as  those  calculated  using  the  kp  k2  combination 
specified  here. 

Pile  drag  coefficients  inferred  from  the  water  channel  tests  were  all 
slightly  less  than  unity.  Thus  a  somewhat  conservative  value  of  1.0  was  chosen 
for  the  simulation.  From  ideal  fluid  flow  theory,  the  added  mass  coefficient 
should  be  2.0.  This  was  increased  to  3.0  to  include  the  additional  effective 
pile  mass  due  to  flooding  of  the  hollow  pile. 
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Table  6.  Parameters  for  full  scale  prototype  modeling. 


Parameter 

OPPILE 

High  Low 

Water  Resonant  Water 

HURPILE 

High  Storm 

Water  Surqe 

COLPILE/RECPILE 

Typical  Barqe 

Orag  coef.  of  boards 

1.28 

Drag  coef.  of  pile 

1.0 

1.0 

1.0 

Added  mass  coeff.  of 

pile 

3 

3 

3 

kl 

577,700  ft- lbs/rad 

577,700  ft-lbs/rad 

577,700  ft-lbs/rad 

k2 

28,900  ft- lbs/rad 

28,900  ft-lbs/rad 

28,900  ft-lbs/rad 

Breakpoint  angle 

10  deg 

10  deg 

10  deg 

Load  weight 

1,744  lbs 

1,744  lbs 

1,024  lbs 

Pile  weight 

2,890  lbs 

2,890  lbs 

2,716  lbs 

Pile  length  (lp) 

40.7  ft 

40.7  ft 

38.25  ft 

Pile  diameter 

1.5  ft 

1.5  ft 

1.5  ft 

Length  to  load  (1^) 

8.84  ft 

8.84  ft 

12.07  ft 

Length  to  boards  (1b) 

37.5  ft 

Area  boards 

36  ft2 

Depth  to  hinge  (d) 

27.5  ft  27.5  ft  7.5  ft 

27.5  ft  36.5  ft 

25  ft 

Total  depth  (dt> 

30  ft  30  ft  10  ft 

30  ft  39  ft 

30  ft 

Wind  velocity 

60  kts  =  101.3  ft/s 

100  kts  =  168.8  ft/s 

Current  velocity 

3  kts  «  5.063  ft/s 

3  kts  =5.063  ft/s 

Wave  height 

5  ft 

6  ft 

Wave  period 

5  s  2.8  s  5  s 

5  s 

Wave  length 

118  ft  41  ft  82  ft 

118  ft  123  ft 

Static  angle 

5.0  deg  5.0  deg  4.2  deg 

4.3  deg  3.8  deg 

Natural  period 

2.83  s 

2.83  s 

2.52  s 

Barge  speed  (Ub) 

10  kts  =  16.88  ft, 

Length  (L) 

180  ft 

Freebard  (fb) 

0 

Bow  angle  (ef) 

25  deg 

Oraft  (dh) 

12  ft 

OPPILE  and  HURPILE  were  used  to  model  operating  and  storm  environments, 
respectively.  OPPILE  computer  simulations  were  conducted  for  high  water/ 
design  wave,  high  water/ resonant  wave  and  low  water/design  wave  conditions. 

For  each  wave  forcing,  the  computer  simulation  was  run  beyond  the  point  for 
which  initial  transients  had  completely  decayed  (after  3  wave  periods)  in 
order  to  evaluate  the  regular,  periodic  response.  Motion  about  the  lower  hinge 
was  considered  since  this  presents  the  largest  projected  area  and  is  therefore 
the  "worst  cast"  for  wave/current  loading. 

COLPILE  and  RECPILE  were  used  to  estimate  the  pile  response  to  a  direct 
barge  hit  at  the  design  speed  of  10  kts.  Motion  about  the  upper  hinge  was 
simulated  here  because  it  represents  the  "worst  cast"  for  barge  clearance  and 
hence  impact  loads.  In  a  representative  simulation,  a  barge  bow  angle  of  25  deg 
and  a  draft  of  12  ft  were  used  since  these  values  are  typical  of  barges  using 
the  Houston  area  ship  channel.  To  determine  the  effect  of  hinge  clearance  on 
impact  loads,  an  additional  series  of  COLPILE  runs  were  made  for  the  possible 
range  of  barge  drafts. 

PREDICTIONS 

The  OPPILE  computed,  steady  state  response  time  series  for  each  of  the  3 
operating  condition  cases  considered  are  shown  in  Fig.  26.  The  motion  due  to 
the  resonant  wave  excitation,  as  expected,  has  the  largest  amplitude  but  does 
not  exceed  the  maximum  allowable  vertical ity  limit  of  10  deg.  At  the  other 
extreme,  the  low  water  response  is  much  reduced  and  less  regular  due  to  weaker 
wave  forcing.  Hinge  moments  necessary  to  maintain  the  less  stringent  of 
the  verticality  criteria  are  seen  to  approach  100,000  ft-lbs. 

The  HURPILE  calculated  steady  state  response  for  the  hurricane  design 
conditions  are  shown  in  Fig.  27.  Predicted  average  angular  position  is  somewhat 
less  and  maximum  angle  amplitude  does  not  greatly  exceed  that  predicted  for 


Time  (s) 


Pile  motion  during  hurrlcan  conditions.  Shown  Is 
the  steady-state  response  over  a  wave  period  as 
computed  by  HURPILE.  Plotted  are  response  for  normal 
high  water  (30  ft.  depth)  conditions  (-)  and  storm 
surge  (39  ft.  depth)  conditions  (...).  The  time 
scale  starts  at  the  beginning  of  a  wave  crest. 


operating  conditions.  This  is  explained  largely  by  the  decrease  in  forcing 
accompanying  the  sacrifice  of  the  daymark  boards. 

Results  calculated  by  COLPILE  for  the  " typical"  barge  collision  are 
given  in  Table  7.  Hinge  moment  and  reaction  force  loads  in  general  seem 
reasonable,  but  the  force  impulse  at  the  bottom  of  the  bow  rake  actually  re¬ 
presents  very  large  impact  forces  due  to  the  short  duration  of  the  impact 
process.  Fig.  28  shows  how  the  lower  bow  rake,  barge  force  impulse  becomes 
even  larger  with  further  decrease  in  clearance.  These  results  suggest  that  it 
will  be  difficult  to  avoid  local  pile  or  hinge  damage  for  clearances  with  respect 
to  the  hinge  less  than  7  ft  unless  special  shock  absorbing  cladding  is  employed. 

Pile  motion  during  recovery  from  the  design  barge  collision,  as  computed 
by  RECPILE,  is  shown  in  Fig.  29.  Returning  to  the  vertical  takes  less  than 
10  s.  The  small  overshoot  and  oscillation  about  the  equilibrium  point  is 
quickly  damped. 

In  general,  it  is  seen  that  CTPS  design  easily  satisfies  the  minimum 
design  criteria  for  worst  case  operating  and  hurricane  conditions.  It  does  not, 
however,  appear  feasible  to  design  CTPS  systems  to  sustain  barge  collisions  at 
very  low  clearance  distances.  With  this  single  qualification,  the  predicted 
performance  has  been  found  to  be  entirely  acceptable. 


Table  7.  Collision  of  "typical  barge"  with  pile.  Position  and  dynamic  loads 

are  computed  by  COLPILE.  Time  is  given  in  s,  8  in  deg,  moments  in  ft- lbs 
forces  in  lbs.  During  impact,  moment  impulse  (ft-lb-s)  and  force  impulse 
(Ib-s)  are  listed. 


Time 

6 

_ 1_ 

fb 

"n 

rh 

RV 

Impact  at  top  of  bow  rake 

0 

0 

68,000 

2,720 

-850 

0 

.25 

9.6 

137,900 

5,470 

96,600 

4,380 

3,060 

.50 

18.7 

124,300 

4,730 

105,200 

4,340 

3,940 

Impact  at  bottom  of  bow  rake 

.65 

25.0 

61,170 

4,260 

1,080 

500 

.75 

28.5 

189,000 

12,840 

110,200 

12,970 

7,050 

1.00 

41.0 

157,300 

9,180 

116,400 

9,490 

9,320 

1.25 

50.0 

137,200 

6,810 

121,000 

6,620 

9,690 

1.50 

56.6 

81 ,500 

3,470 

124,300 

3,580 

7,990 

1.75 

61.5 

62,700 

2,310 

126,800 

2,270 

7,080 

2.00 

65.2 

57,600 

1,870 

128,700 

1,600 

6,530 

2.25 

68.1 

57,400 

1,650 

130,150 

1,220 

6,160 

Contact  with  barge 

bottom 

2.5-12.8 

70.1 

70,700 

1,910 

131,200 

190 

5,640 

68 


Hinge  Clearance  (ft) 


Fig.  28.  Barge  force  Impulse  as  a  function  of  barge  draft 
(or  hinge  clearance).  Total  water  depth  Is  30  ft 
and  hinge  height  Is  5  ft. 


Pile  recovery  as  computed  by  RECPILE.  Shown  are 
pile  recovery  from  a  "typical  barge"  collision  (  — ) 
and  recovery  from  a  90  deg.  knockdown  in  an  adverse 
current  of  3  knots  (...). 


IX.  LOGISTICAL  ANO  COST  FACTORS 

SOIL  MECHANICS 

The  ability  of  the  subbottom  sediment  to  withstand  the  loading  due  to  the 
impact  process  is  crucial  to  the  success  of  the  CTPS  system.  The  base  must  remain 
fixed  in  order  for  the  design  criteria  to  be  met.  What  follows  is  a  first  order 
analysis  of  the  soil  mechanics  problem  which  will  provide  sufficient  information 
for  practical  design  decisions. 

The  problem  is  that  of  a  circular  pile  driven  into  the  ocean  bottom  subjected 
to  a  laterally  applied  impact  force.  The  conventional  approach  to  laterally 
loaded  piles,  however,  is  a  static  analysis,  and  this  method  is  used  here.  It  is 
assumed  that  the  sediment  resistance  varies  linearly  with  depth  from  a  value  of 
zero  at  the  sediment  surface  to  a  maximum  value  at  the  bottom  of  the  pile.  Using 
a  graphical  solution  presented  in  Fig.  72  of  Paules  and  Davis  (1980),  and  assuming 
a  typical  value  of  bouyant  density  (total  sediment  density  minus  the  water  density) 
at  60  Ib/ft  and  a  pile  diameter  of  1.5  ft,  results  for  horizontal  force  R^  on  a 
30  ft  and  a  40  ft  embedded  pile  section  are  8,300  lbs  and  11,900  lbs,  respectively. 

These  forces  are  significantly  smaller  than  an  estimated  hinge  reaction  force 
during  impact  of  27,000  lbs  calculated  using  COLPILE  and  collision  experiment 
results.  Paules  and  Davis  indicate  that  there  are  various  design  modifications 
for  correcting  this  problem  and  effectively  strengthening  the  pile  system  with 
respect  to  laterally  applied  loads.  The  basic  concept  inherent  in  all  the 
schemes  is  to  increase  the  soil  bearing  surface. 

The  most  desirable  solution  is  to  add  wings  to  the  upper  portion  of  the 

driven  pile.  The  wings  can  be  analyzed  as  retaining  walls  in  a  passive  earth 

pressure  case  under  undrained  conditions.  Employing  the  method  presented  by  Lambe 

3  2 

and  Whitman  (1969)  and  using  60  Ib/ft  for  bouyant  density  and  500  lb/ft  for 
ultimate  shear,  a  lateral  force  of  over  34,500  lbs  can  be  sustained  by  a  3  ft 
wide  by  5  ft.  high  rectangular  wing.  Four  wings  of  this  size  could  be  added  to 


the  upper  portion  of  the  driven  pile  to  achieve  more  than  ample  resistance  to  the 
applied  lateral  load. 

This  analysis  indicates  that  a  stand  alone  single  pile  cannot  withstand  the 
laterally  applied  impact  load.  The  addition  of  moderately  sized  wings,  on  the 
other  hand,  appears  to  be  a  technically  feasible  solution.  More  analysis  con¬ 
cerning  the  site  specific  soil  profile  is  necessary,  however,  and  should  be 
completed  for  each  installation. 

SURVEY  OF  HOUSTON  AREA  NAVIGATION  AIDS 

A  survey  was  taken  of  the  navigation  aids  in  the  Houston  area  to  evaluate 
the  extent  to  which  existing  markers  are  replacable  by  CTPS  systems.  The 
geographical  region  considered  is  the  Houston  Ship  Channel  starting  at  Point 
Bolivar  and  extending  across  Galveston  Bay  to  Morgan  Point.  According  to  the 
Coast  Guard  "Light  List"  (1984),  this  area  includes  a  total  of  74  in-water, 
government  navigation  aids. 

To  determine  which  aids  may  be  exchanged,  criteria  based  on  water  depth 
were  developed.  It  was  assumed  that  a  CTPS  can  sustain  a  collision  at  maximum 
speed  by  a  fully  loaded,  typical  barge  (draft  =  12  ft)  if  the  fraction  of  water 
depth  required  for  clearance  is  the  same  as  in  the  design  case.  Since  5  ft  for 
hinge  height  plus  7  ft  to  reduce  impact  forces  provides  a  generous  clearance  in 
30  ft  of  water,  the  minimum  depth  can  be  evaluated  from 

(>5)/30  =  (16) 

Thus  the  minimum  depth  is  (d*)„,.  =  20  ft.  A  maximum  depth  of  40  ft  was  chosen 

t  min 

as  an  upper  limit  since  the  CTPS  design  can  easily  be  scaled  upwards  by  1/3. 

The  number  of  in-water  government  aids,  including  buoys,  dolphins  and  piles,  in 
the  depth  range  20  to  40  ft  is  14. 

This  criteria  may  be  relaxed  if  It  is  assumed  that  the  aid  Is  replaceable 
for  depths  In  which  the  typical  barge  must  simply  clear  the  hinge  itself.  Since 
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the  absolute  minimum  mud  clearance  for  a  30  ft  water  depth  is  6.5  ft.,  the 
less  restrictive  minimum  depth  criterion  becomes 

6.5/30  =  [<dt),ln  -  12Kdt)„in  (17) 

yielding  (dt)M1-n  =  15  ft.  This  criterion  should  be  applicable  if  the  barge  does 
not  hit  directly  at  full  speed  or  if  the  barge  is  not  fully  loaded.  The  number  of 
in-water  government  aids  (excluding  skeleton  towers)  in  the  depth  range  from  15 
to  19  ft  is  22. 

In  summary,  14  out  of  74  or  19%  are  replaceable  using  a  criteria  based  on 
the  maximum  design  collision  condition  with  a  typical  barge.  Nearly  half  (14  + 

22  j=  36  out  of  74),  however,  are  replaceable  assuming  that  only  minimum  clearance 
for  a  fully  loaded  barge  is  necessary. 

INSTALLATION  PROCEDURES 

Installation  of  the  CTPS  should  not  impose  serious  departures  from  present 
practices  of  installation  and  does  not  require  the  use  of  divers  or  sophisticated 
underwater  equipment  as  required  by  the  logistical  design  criteria  given  in 
Table  1.  Use  of  an  18  inch  steel  pipe  as  piling  is  familiar  to  the  Coast  Guard 
since  it  is  presently  in  use  for  ordinary  pile  structures.  The  hinge  component 
can  be  pre-assembled  on  shore  to  reduce  crew  effort  and  vessel  time  during  the 
actual  deployment. 

The  total  length  of  the  pile  system  is  44  ft  plus  at  least  30  ft  below 
the  hinge.  The  embedment  length  of  30  ft  normally  needed  to  obtain  fixity  in 
soft  clay  may,  however,  have  to  be  increased  because  of  the  hinge  impact  loads 
transmitted  to  the  base.  Since  the  total  length  exceeds  the  50  ft  maximum 
length  stipulated  by  requirement  (b)  in  Table  1  (due  to  limited  crane  height), 
the  pile  must  be  installed  in  sections. 

A  scheme  for  dividing  the  CTPS  into  3  parts  is  shown  in  Fig.  30.  The 
above  water  part  is  attached  via  a  flange  coupling  to  the  main  pile  which  in- 


eludes  the  spring,  hinge  and  base.  The  bottom  of  the  base  consists  of  a  tapered 
socket  which  fits  over  the  top  of  the  lowest  section  embedded  in  the  soil. 

Gussets  are  used  to  reinforce  the  base/socket  joint  and  to  increase  soil  bearing 
area.  The  enlarged  soil  contact  surface  aids  the  embedded  pile  section  in  re* 
si  sting  base  movement  during  collisions. 

The  installation  procedure  begins  by  driving  the  lower  section  until  its 
top  is  just  above  the  water  surface.  (Thus  its  length  is  governed  by  water 
depth  and  may  be  longer  than  that  necessary  to  obtain  fixity).  Next  the  main 
pile  member,  complete  with  a  full-length  driving  collar,  is  set  in  place  and 
the  socket  joint  secured  by  bolts.  The  fastening  need  not  be  extensive  since 
the  socket  itself  is  designed  to  transfer  the  bending  moment  and  compressive 
axial  and  shear  forces  to  the  lower  section.  The  bolts  are  mainly  a  precaution 
against  minor,  inadvertent  lifting  which  might  occur  during  installation. 

Using  the  collar  to  transfer  the  hammer  blows  directly  to  the  base  plate, 
the  system  is  driven  so  that  the  base  is  level  with  the  mudline.  After  driving, 
the  (reuseable)  collar  is  removed.  Lastly  the  top  section  is  attached,  stay 
tension  is  adjusted  and  the  navigation  aid  is  built.  For  shallow  areas,  of 
course,  the  shorter  CTPS  sections  could  be  joined  on  deck. 

When  removing  the  CTPS,  the  structure  will  separate  at  the  base  socket/lower 
pile  joint.  Thus  the  major  portion  of  the  system  is  recovered  intact,  and  the 
remaining  section  is  below  the  mudline.  The  bolting  used  to  secure  this  joint 
during  installation  will,  therefore,  be  done  so  that  it  serves  as  a  weak  point 
with  respect  to  direct  upward  pull. 

COST  ESTIMATE 

Cost  estimates  (materials  plus  fabrication)  were  made  for  the  3  major  design 
components  -  the  piping,  the  hinge/stay  system  and  the  "spring"  element.  Sufficient 
piping  for  both  the  section  above  the  hinge  and  the  embedded  portion  is 
estimated  to  cost  $4,000.  Hinge/stay  system  construction  should  total  approxl- 
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mately  $2,000,  while  "spring"  costs  may  vary  from  $800  to  $2,900  depending  on  the 
type  of  “spring"  used. 

Complete  system  costs  will  therefore  be  between  $6,800  and  $8,900.  Using 
information  supplied  by  Miller  (1982),  installation  costs  are  approximately 
$4,000,  so  total  costs  range  from  $10,800  to  $12,800.  In  comparison  the  total 
expense  (pile  plus  installation)  of  a  simple  12  inch  wooden  pile  is  about  $4,500. 
Thus  the  CTPS  approach  would  become  cost-effective  after  3  collisions. 


X.  DISCUSSION 


CONCLUSIONS 

Bench  testing  of  the  CTPS  design  developed  in  this  study,  performed  out  of 
water  using  the  physical  scale  model,  showed  first  of  all  that  the  hinge  possesses 
the  necessary  kinematic  characteristics.  Full  90  deg  from  the  vertical  articulation 
is  possible  for  all  horizontal  angles.  Secondly  the  desired  piecewise  linear  hinge 
moment  as  a  function  of  angle  behavior  (M^  =  (r))  was  observed.  The  scale  model 

experiments  indicated  that  the  design  initial  stiffness  coefficient  k^  =  11.4 
ft-lb/rad  (577,700  ft-lb/rad  full-scale)  can  be  achieved,  the  actual  breakpoint 
angle  occurs  approximately  at  the  design  value  of  10  deg,  and  normally  k2  s  1/20 
k-|.  It  was  noted  that,  depending  on  the  prestressed  stay  tension  balance  and 
internal  friction,  k-j  can  exceed  the  design  value  though  large  angle  hinge 
moments  are  not  greatly  affected.  This  increase  in  initial  stiffness,  however, 
enhances  the  CTPS  operational  vertical ity  performance  and  is  not  viewed  as 
detrimental. 

Water  channel  experiments  conducted  at  the  Coast  Guard  Academy  CWT  facility 
showed  that,  at  the  model  scale  (1/15)  the  design  easily  meets  the  vertical ity 
requirement  at  the  maximum  current  specified  in  the  design  criteria.  Lateral 
movement  due  to  vortex  shedding  was  observed  to  be  negligible.  Measured  effective 
drag  coefficients  for  the  pile  in  water  were  found  to  be  less  than  1  so  that  use 
of  Cw  «  1  for  computer  modeling  purposes  may  be  viewed  as  being  conservative. 

Collision  experiments,  using  the  physical  scale  model  in  the  UNH  indoor  pool, 
showed  that  the  CTPS  design  can  sustain  impact  at  the  required  (scaled)  barge 
speed  and  recover  to  an  upright  position.  Impact  forces  were  measured,  and 
force  impulses  were  compared  with  predictions  made  using  the  computer  models. 
Agreement  was  within  23%  indicating  the  computer  simulations  are  sufficiently 
accurate  for  design  evaluation  purposes. 
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Computer  simulations  were  made  for  operational,  hurricane  and  collision 
conditions  using  imput  dimensions  from  the  design  itself,  performance  parameters 
found  suitable  from  the  physical  model  experiments,  and  environmental  specifications 
required  by  the  design  criteria.  The  computer  program  OPPILE,  developed  to 
determine  pile  performance  during  operating  conditions,  showed  that  the  design 
meets  the  maximum  allowable  verticality  limit.  The  computer  model  developed  for 
calculating  pile  motion  during  hurricane  conditions,  HURPILE,  indicates  that, 
since  the  daymark  boards  are  sacrificed,  pile  movements  and  hinge  moments  are  not 
large,  and  no  damage  should  be  expected.  The  programs  developed  for  modeling 
collision  processes,  COLPILE  and  RECPILE,  demonstrate  that,  for  collisions 
at  the  design  speed  by  typical  barges  used  in  the  Houston  area,  damage  to  the  pile 
itself  should  not  occur  and  pile  recovery  is  expected  to  be  prompt.  COLPILE  does 
predict,  on  the  other  hand,  that  impact  of  the  barge's  bottom  of  the  bow  rake 
with  the  pile  has  the  greatest  potential  for  damage.  This  problem  increases 
in  severity  with  decrease  in  clearance.  Because  of  the  need  for  sufficient 
clearance  and  the  height  required  for  the  hinge  itself,  the  very  severe  mud 
clearance  of  3  ft  specified  in  the  design  criteria  cannot  be  met.  As  noted 
above,  however,  there  should  be  no  difficulty  with  collisions  by  barge  traffic 
of  average  draft. 

Logistical  requirements  for  the  design  were  satisfied  since  installation 
procedures  do  not  differ  radically  from  existing  practice.  A  method  of  in¬ 
stalling  the  pile  in  sections  was  found  to  be  suitable  for  meeting  the  pile 
length  limitation  of  50  ft.  Costs  for  the  design  were  estimated  to  fall  between 
the  design  goal  of  $5,000  and  the  maximum  allowable  amount  of  $10,000.  At  this 
level  of  expense,  a  deployed  CTPS  will  become  cost-effective  after  three 
collisions. 

At  the  present  time,  a  1/4  scale  model  is  being  built  for  field  testing  using 
Sea  Grant  funding.  It  is  expected  that  this  experiment  will  shed  additional  light 


on  many  of  the  practical  engineering  aspects  such  as  material  selection,  system 
fabrication  and  installation. 

RECOMMENDATIONS  ■ 

The  pile  system  design  developed  in  this  study  appears  to  be  a  feasible 
solution  to  the  problem  of  pile  damage  in  the  nation's  extensive  high  traffic, 

i 

shallow  waterways.  Physical  scale  models  and  computer  simulations  indicate  that 
the  design  is  technically  sound  and  preliminary  cost  estimates  have  shown  that  the 
concept  is  financially  beneficial  where  the  potential  for  collision  is  great. 

i 

We  therefore  recommend  that  the  next  step  in  development  be  taken  in  which  a 
full-scale  prototype  is  built  and  field  tested.  Construction  plans  for  prototype 
fabrication  should  be  based  principally  on  the  design  presented  in  this  report. 

i 

Departures  to  conform  to  a  particular  installation  location,  of  course,  may  be 
made.  It  is  recommended,  however,  that  design  parameters  be  re-established  using 
the  design  guidelines  and  computer  models  as  well  as  the  physical  scale  model 

i 

data  presented  here. 

In  summary,  the  results  obtained  in  this  study  are  very  encouraging.  Con¬ 
tinued  development  and  subsequent  implementation  of  working  systems  is  recommended 
and  promises  eventually  to  save  much  effort  and  expense  on  the  part  of  the  Coast 
Guard. 
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VII.  APPENDIX  A:  COMPUTER  PROGRAM  LISTINGS 


STATUS  OF  PROGRAMS 

In  this  appendix  program  listings  for  the  computer  models  described  in 
Section  II  are  provided.  The  programs  are  written  in  Applesoft  BASIC  and  were 
originally  used  on  an  Apple  lie  microcomputer  with  printer.  The  type  and  method 
of  input,  pre-set  parameter  values  (changeable  by  editing),  and  output  options 
are  evident  from  the  listings. 

The  programs  are  presented  here  as  they  exist  at  the  time  of  the  writing  of 
this  report  (November  1984)  though  they  are  presently  being  upgraded.  Anticipated 
changes  include  improving  program  efficiency  and  making  the  programs  clearer 
to  potential  users.  Though  the  programs  are  believed  fundamentally  correct, 
further  verification  studies  are  also  planned.  The  authors  would  appreciate 
being  informed  of  any  ''bugs"  or  errors  found  by  others.  Due  to  the  transitional 
status  of  these  programs,  the  authors  cannot  be  responsible  for  their  misuse. 

No  warranty  is  expressed  or  implied. 


PILESTIFF 


100 

PRINT  “INITIAL  STIFFNESS  FOR  STATIC 

CONDITIONS 

103 

PRINT 

110 

INPUT  “INPUT 

THETA  (DEC)  “-.THETA 

US 

PRINT 

120 

INPUT  "INPUT 

LENGTH  TO  LOAD  (FT)  “; 

LM 

123 

PRINT 

130 

INPUT  “INPUT 

WOT  OF  LOAD  (LSS)  ":WL 

133 

PRINT 

140 

INPUT  “INPUT 

PILE  LENGTH  (FT)  “;LP 

143 

PRINT 

ISO 

INPUT  "INPUT 

PILE  VCT  (LSS)  “:VP 

133 

PRINT 

1(0 

INPUT  “INPUT 

DEPTH  TO  HINGE  ":D 

US 

PRINT 

170 

INPUT  “INPUT 

PILE  01 AN  (FT)  *:DP 

173 

PRINT 

ISO 

INPUT  "INPUT 

VINO  VEL  (FT/S)  ":UA 

103 

PRINT 

190 

INPUT  “INPUT 

LENGTH  TO  BOARDS  (FT) 

“ ;  LS 

193 

PRINT 

200 

INPUT  “INPUT 

AREA  OF  BOARDS  (FTa2) 

“ ;  AS 

203 

PRINT 

210 

INPUT  “INPUT 

CUR  VEL  (FT/S)  “:UC 

213 

PRINT 

230  aw  .  (4.0  /  32.2 
240  CA  •  1.0 
230  CV  a  1  .0 
2(0  CS  a  1.20 

270  XI  a  LN  *  VL  ♦  .3  •  LP  »  VP  ♦  (.23  «  HA  •  CA  *  OP  •  <L?  a  2  -  0  a 
2)  •  <UA  a  2)  ♦  .3  *  BA  •  C8  *  AB  »  LS  *  (UA  A  2)  ♦  .  23  •  (0  a  2)  «  RV 
*  CV  *  OP  *  <UC  A  2>>  /  (THETA  *  3.1414  /  100) 

200  PRINT  "INITIAL  STIFFNESS  XI  (FT-LBS/RAD)  a  ";X1 
290  END 


-I*. .MV1.1  --'j • . • . - 

PILEFREQ 

P.-V 

100 

PRINT  “7IL2  NATURAL  FREQUENCY" 

tos 

PRINT 

• 

110 

INPUT  "INPUT  INIT  STIFF  XI  (FT-LBS/RAD)  ";X1 

/* 

us 

PRINT 

120 

INPUT  "INPUT  LENGTH  TO  LOAD  (FT)  ";LH 

•v  >; 

•  V,‘ 

123 

PRINT 

*  a"  a 

130 

INPUT  "INPUT  VCT  OF  LOAO  (LSS>  ":WL 

133 

PRINT 

• 

no 

INPUT  "INPUT  LENGTH  OF  FILS  (FT)  ";LP 

143 

PRINT 

130 

INPUT  "INPUT  VGT  OF  FILS  (LBS)  ";VP 

133 

PRINT 

140 

INPUT  "INPUT  SIAN  OF  FILE  (FT)  «:OP 

»  - 

US 

PRINT 

i 

170 

INPUT  "INPUT  DEPTH  TO  HINGE  (FT)  ":D 

173 

PRINT 

100 

RV  •  44.0  /  32.2 

190 

CM  •  2.0 

200 

HP  >  VP  /  32.2 

fc  .-*=■-  A 

210 

ML  -  VL  /  32.2 

1 

220 

IK  a  HP  *  (LP  a  2)  /  3  ♦  (Lrt  a  2)  *  NL  ♦  (3.1414  /  12)  *  HV  ■  CH  * 

V*\ 

I  (OF  a  2)  «  (0  a  3) 

230 

NF  «  (<K1  -  LH  *  VL  -  .3  *LP  *VP)  /  IH>  a  .3 

-  ‘.V 

240 

PRINT  "NATURAL  FREQUENCY  (RAO/S)  a  «;NF 

•  *'*.’•  . 

243 

PRINT 

230 

PRINT  "NATURAL  PERIOD  <S)  a  ":2  *  3.X4I4  /  NP 

^  . - >4 

240 

END 

•**.  •* 

. 

»"  / 

■*.  •*.  •* 

v'v 

» 

'  .  •  ,  « 

1 _ 

.*’»/*•* 

>  .*•  .*• 

A- 3 

1 _ 

>  1*1 

WAVELENGTH 


;o 

PRINT 

“WAVELENGTH  CALCULATION” 

20 

PRINT 

“FOR  A  SECULAR  WAVE” 

23 

PRINT 

10 

INPUT 

“INPUT  WAVS  PERIOD  <S)  “:?ER 

33 

PRINT 

40 

INPUT 

“INPUT  DEPTH  <FT)  ":D 

43 

PRINT 

30 

PI  •  3. 

1414 

40  G  •  32 . 2 
70  i  •  .01 

lo  to  -  c  *  <pes  a  ])  /  (i  *  ;n 

*0  L  >  to 

100  F  a  <L  /  LO)  -  (  El?  (2  •  PI  *  0  /  L)  -  XXP  (  -  2  •  FI  *  D  /  L)> 

/  <  EXP  <2  *  PI  »  0  /  t)  ♦  EXP  <  -  2  *  PI  *  0  /  L>> 

120  IF  ASS  (F)  <  X  THEN  GOTO  140 

130  OF  •  (1  /  LO)  ♦  ((  EXP  (2  *  PI  *  0  /  L)  +  EXP  (  -  2  *  PI  *  0  /  L) 

)  a  <  -  2) )  *  0  *  PI  «  0  /  <L  a  2) 

HO  L  »  L  -  F  /  OF 
ISO  GOTO  100  . 

140  PRINT  “WAVELENGTH  •  “:L 


OPPILE 


1  PRINT  “PILE  MOTION  SOS  QP5HATI ONAL  CONDITIONS" 

5  PUNT 

10  INPUT  “INPUT  STEPS I ZE  (S)  " :H 
IS  PRINT 

20  INPUT  “INPUT  MAX  T  (S)  “ ;MAX 
25  PRINT 

30  L«  •  INT  (MAE  /  H> 

40  INPUT  “INPUT  1ST  STIFFNESS  CONSTANT  El  (FT-L8S/RAD)  “;E1 
45  PRINT 

SO  INPUT  -INPUT  2ND  STIFFNESS  CONSTANT  X2  (FT-LBS/RAD)  “:K2 
55  PRINT 

40  INPUT  -INPUT  BREAKPOINT  ANCLE  (DEC)  “  :8A 

41  B THETA  a  BA  *  3.1414  /  ISO 
45  PRINT 

70  INPUT  “INPUT  LOAD  VGT  UBS)  “;VL 

75  PRINT 

10  ML  a  Vt  /  32.2 

♦0  INPUT  “INPUT  PILE  VCT  (LBS)  “;VP 


75 

PRINT 

100 

MP  a  VP  /  32.2 

110 

INPUT  “INPUT 

PILE  LENGTH  (FT)  “UP 

115 

PRINT 

120 

INPUT  “INPUT 

PILE  DIAMETER  (FT)  “;DP 

125 

PRINT 

130 

INPUT  “INPUT 

LENGTH  TO  LOAD  (FT)  “;LN 

133 

PRINT 

140 

INPUT  “INPUT 

LENGTH  TO  BOARDS  (FT)  “;LB 

145 

PRINT 

130 

INPUT  “INPUT 

DEPTH  TO  HINGE  (FT)  “;D 

135 

PRINT 

140 

INPUT  “INPUT 

TOTAL  DEPTH  (FT)  ■ :DT 

143 

PRINT 

170 

INPUT  “INPUT 

VINO  VSL  (FT/S)  “jUA 

173 

PRINT 

100 

INPUT  “INPUT 

CUR  VEL  (FT/S)  “:UC 

IIS 

'PRINT 

170 

INPUT  “INPUT 

VAVE  HGT  (FT)  “;HV 

173 

PRINT 

200 

INPUT  “INPUT 

VAVX  PERIOD  (S)  “ s PER 

203 

PRINT 

210 

SICMA  a  2  *  3 

.1414  /  PER 

220 

INPUT  “INPUT 

VAVE  LENGTH  (FT)  “i LAMBDA 

223  PRINT 

230  XA  a  2  •  3.1414  /  LAMBDA 

240  RA  a  .077  /  32.2 

250  RV  a  44.0  /  32.2 

240  AB  a  34 

270  CA  a  1.0 

210  CV>  t.l 

270  CB  a  1.21 

300  CM  a  2.0 

302  Cl  a  .5  •  RV  «  CV  *  DP 

304  C2  a  .5  *  KV  •  SICHA  /  (  SEP  (XA  *  DT)  -  EXP  (  -  XA  *  DT)) 
304  C3  a  -  (3.1414  /  0)  *  RV  •  CM  »  (DP  a  2)  *  (SICMA  a  2)  •  KV 
XP  (XA  ■  DT)  -  EXP  (  -  XA  •  0T>) 


no  IK  «  <MP  «  L?  A  :j  /  J  +  ML  *  LM  a  2  *  CM  *  RV  *  (3.14U  /  12)  *  ‘ 

0?  A  ;>  *  <0  A  3) 

320  SIM  T(L*  *  2.3) 

330  SIM  THETA (  L%  *  2.3) 

340  SIM  OMEGA ( L%  *  2.3) 

3S0  SIM  EC  3) 

340  SIM  M<  3) 

370  SIM  MHO) 

300  SIM  MC<  3) 

390  DIM  MV<  3) 

400  SIM  MCO) 

410  1%  •  0 
420  T( 0 . 0)  •  0 

430  THETA(  0.0)  •  (.23  *  (LP  a  2  -  0  a  2)  •  HA  •  CA  *  DP  '  «A  a  J  ♦  .3 

*  LB  *  HA  *  CS  *  AB  *  UA  a  2  ♦  .  23  *  D  *  RV  *  CV  *  DP  *  0  »  UC  a  2 )  / 

(HI  -  (LM  *  VL  ♦  .3  *  LP  *  VP)) 

440  PRINT  "STATIC  EOUIL  ANCLE  ■  " ;THETA(0 . 0 )  •  110  t  3.1414 
443  PRINT 

447  PRINT  "TIME  <S>  a  ";T(I%.3);"  THETA  (DEC)  a  “  ;THETA(  1%.  0)  *  100 
/  3.1414 

430  T(I*  ♦  1.0)  a  T( I*. 0)  *  X 

440  II  7(1%  >  1.0)  >  MAX  THEN  GOTO  70S 

470  T(I%.1)  a  T(I*.0)  ♦  .3  *  H 

400  T( 1%. 2)  a  T( I%.l) 

490  T(I*.3)  ■  T( I*. 0)  ♦  X 
300  J*  •  0 

310  K(J«)  «  OMEGA ( 1%. J%) 

320  IF  THETA ( 1%.  J%)  )  (STHETA)  THEN  GOTO  340 
330  IF  THETA(I*.J%)  ♦  1THSTA  <  0  THEN  GOTO  380 
340  MK(J«)  a  XI  *  THETA(I«.J«) 

330  GOTO  390 

34  0  HH(J»)  a  XI  •  HTHETA  ♦  X2  *  (THETAC 1%. J%>  -  HTHETA) 

370  GOTO  390 

300  MH(J%)  a  -  XI  •  HTHETA  ♦  X2  »  (THETA(  I*.  J%>  ♦  HTHETA) 

590  MC(J%)  a  (LM  •  VL  ♦  .3  «  LP  *  VP)  •  SIN  (THETA(  I*.  J%) ) 

400  MV(J»)  .  .23  •  ((LP  *  COS  (THETA(  II. J%)  ) )  a  2  -  D  a  2)  •  HA  »  CA 

•  DP  •  UA  a  2  ♦  Li  *  COS  (THETA(I*.J%>>  •  . 3  •  HA  •  CB  *  AS  •  UA  a  2 

410  DEF  FN  WHY)  *  (D  ♦  Y)  *  Cl  *  ((C2  *  (  EXP  (XA  •  DT  ♦  KA  *  Y)  ♦ 

EXP  (  -  KA  *  OT  -  XA  *  Y))  *  COS  (SIGMA  *  T(I*.J%>>  ♦  UC  -  <D  ♦  Y)  » 
OMEGA  (I*.  J%) )  a  2)  •  SON  (C2  *  (  EXP  (XA  *  DT  ♦  XA  *  Y)  ♦  EXP  (  -  X 
A  *  DT  -  KA  «  Y))  *  COS  (SIGMA  *  T(I*.J%>>  ♦  DC  -  (D  ♦  Y>  *  OMECA(I%. 
J*)) 

412  DEF  FN  HM( Y)  a  (O  ♦  Y)  *  C3  *  (  EXP  (XA  *  DT  ♦  XA  *  Y)  ♦  EXP  ( 

-  KA  *  DT  -  XA  »  Y))  *  SIN  (SIGMA  •  T(X%.J*)> 

420  DLTA  a  0  /  10 

430  SUM  a  0 

440  FOR  N  a  1  TO  9 

430  Y  a  .Da  (DLTA  *  N) 

440  SUM  a  SUM  ♦  FN  HD(Y)  ♦  FN  MM(Y) 

470  NEXT  N 

400  MC(J%>  a  DLTA  *  (SUM  ♦  .3  *  (  FN  MD(  -  D>  ♦  FN  MM(  -  D)  ♦  FN  MD( 
0)  ♦  FN  MM(fl))) 

490  M(J%)  •  (  -  MH( J%)  ♦  MG(J%)  ♦  MV(J%>  ♦  MC(J*))  /  IH 
700  IF  J%  a  3  THEN  GOTO  750 

710  THETA ( 1%.  J%  *  1)  »  THXTACIR.O)  *  (T(I*.J*  ♦  1>  -  T(I*.0)>  •  XCJ%> 
720  OMEGA  ( 1%.  ♦  l)  a  OMEGA  ( 1%.  0)  ♦  (T(I%.J%  ♦  l>  -  T(I%.0)>  *  H(J%> 

730  J»  a  J%  ♦  l 
740  GOTO  310 

730  THETACI*  ♦  1.0)  a  THSTA(I%.0)  ♦  (H  /  4)  *  (X(0>  ♦  2  •  K(l)  ♦  2  •  X 
(2)  a  K(3) ) 


740  OMEGA < !%  ♦  1.0)  *  OMEGA ( 1%. 3)  ♦  CK  /  4)  *  <M(0>  ♦  I  «  M(l>  ♦  2  *  M 
( 2 )  ♦  M(3 ) i 
770  1%  m  1%  *  i 

775  PRINT  "TIME  (3)  ■  •  :T<  I*.  0)  ;  •  THETA  (DEC)  »  “ : THETA <  1%.  0>  «  100 

/  3.1414 

710  GOTO  450 

70S  PRINT 

770  INPUT  “INPUT  INTERVAL  (•  OP  TIME  STEPS)  FOR  OUTPUT  " ;  ITRVL 
000  PRt  i 

iio  i*  «  o 

120  IP  THETA! 1%. 0)  >  < BTHETA )  THEN  GOTO  340 
t30  IP  THETA!  I*. 0)  ♦  BTHETA  <  0  THEN  GOTO  100 
140  MH! 0)  -  XI  *  THETA! I*. 0) 

•50  GOTO  170 

140  MH! 0)  «  XI  *  BTHETA  ♦  X2  *  (THETA!  I*.  0)  -  BTHETA) 

•70  GOTO  170 

•00  MH!  0 )  a  -  XI  *  BTHETA  ♦  X2  »  (THETA!  I*.  0)  ♦  BTHETA) 

•  70  PRINT  “T  !S)  .  "  :T!  1%.  0 ) ; "  THETA  !OEG)  a  “  ;  THETA!  I%.0)  *  180  /  3. 
1414;“  HINGE  MOMENT  !PT-LSS)  a  ";MH!0> 

•75  PRINT 

700  I*  a  1«  ♦  ITRVL 

710  IP  I«  )  L%  THEN  GOTO  730 

720  GOTO  (20 

730  PR»  9 

740  END 


HURPILE 


I  PRINT  "PILE  MOTION  FOR  HURRICANE  CONDITIONS" 
S  PRINT 


40  INPUT  -INPUT  1ST  STIFFNESS  CONSTANT  XI  (FT-tSS/RAO)  -;XI 


SO  INPUT  -INPUT  2ND  STIFFNESS  CONSTANT  X2  (FT- LBS /RAD)  -;K2 


10 

INPUT 

13 

PRINT 

20 

INPUT 

23 

PRINT 

30 

L%  a  I 

40 

INPUT 

43 

PRINT 

SO 

INPUT 

S3 

PRINT 

40 

INPUT 

41 

BTHETA 

43 

PRINT 

70 

INPUT 

73 

PRINT 

10 

ML  a  W1 

70 

INPUT 

73 

PRINT 

100  HP  a  1 

40  INPUT  -INPUT  3REAXPOINT  ANCLE  (DEC)  - : BA 


70  INPUT  -INPUT  PILE  WCT  (LBS)  -;VP 


110  INPUT  "INPUT  PILE  LENGTH  (FT)  ":LP 
US  PRINT 

120  INPUT  -INPUT  PILE  DIAMETER  (FT)  ":DP 
1 25  PRINT 

130  INPUT  -INPUT  LENGTH  TO  LOAD  (FT)  -jLM 
1 35  PRINT 

130  INPUT  -INPUT  DEPTH  TO  HINGE  (FT)  -;D 
133  PRINT 

140  INPUT  -INPUT  TOTAL  DEPTH  (FT)  *:DT 
US  PRINT 

170  INPUT  -INPUT  WIND  VEL  (FT/S)  ";UA 
173  PRINT 

100  INPUT  -INPUT  CUR  VEL  (FT/S)  ";UC 
103  PRINT 

170  INPUT  "INPUT  WAVE  HOT  (FT)  -;HW 
173  PRINT 

200  INPUT  -INPUT  WAVE  PERIOD  (S)  ”;PER 
203  PRINT 

210  SIGMA  >  2  *  3.1414  /  PER 

220  INPUT  -INPUT  WAVE  LENGTH  (FT)  -'.LAMBDA 

223  PRINT 

230  XA  a  2  *  3.1414  /  LAMBDA 

240  RA  a  .077  /  32.2 

230  RW  a  44.0  /  32.2 

240  AB  •  34 

270  CA  •  1.0 

200  CV  a  1  .0 

300  CM  a  2.0 

310  IM  a  (MP  *  LP  a  2)  /  3  a  ML  *  LM  a  2 

320  DIM  T(L%  a  2.3) 

330  DIM  THETA ( L%  a  2.3) 

340  DIM  OMEGA (L%  a  2,3) 

330  DIM  X( 3) 

340  DIM  N( 3) 

370  DIM  MH(3> 


v*'. V* Vv vV* V\. % . ** V “*,* v  v  **\ 


%  s 


380  DIN  MG(3) 

390  OIK  MV(3) 

400  OIK  KC(3) 

410  I»  •  0 
420  T(O.O)  •  0 

430  THETA!  0.0)  •  (.23  *  <1?  a  2  -  D  a  2)  «  RA  *CA  *  DP  '  W  *  1  ♦  .5 

•  LB  »  RA  •  CB  *  AS  •  UA  a  2  ♦  .23  *  0  *  RV  *  CV  *  OP  *  0  *  UC  a  2)  / 

(Xi  -  (iM  •  va  ♦  .3  •  it  *  vp>) 

440  PRIHT  “STATIC  SOUIl  ANGIE  .  “ ; THETA (0.0)  •  180  /  3.1414 
443  PRINT 

447  PRINT  “TINS  (S)  •  “:T(I%.0);“  THETA  (DEC)  >  “ :TH£TA( 1%. 0)  «  180 
/  3.1414 

430  T(I%  *  1.0)  •  T(I«.0)  ♦  H 

440  IP  T(I%  >1.0)  )  MAX  THEN  GOTO  783 

470  T(I%.1)  •  T(I».0)  ♦  .3  •  H 

480  T(I%.2)  -  T(I*.t) 

490  T(I».3>  •  T(I«.0)  ♦  H 

300  J%  «  0 

310  X(J%)  «  OMEGA ( 

320  IP  THETA ( 1%, J%)  )  (BTHETA)  THEN  GOTO  340 
330  IF  THETA ( 1%. J%)  ♦  STHETA  (  0  THEN  GOTO  380 
340  MH(J%)  »  XI  *  THETA ( 1%. J%) 

330  GOTO  390 

340  MH(J%)  >  XI  •  STHETA  ♦  X2  *  ( THETA ( 1%. J%)  -  BTHETA) 

370  GOTO  390 

380  MH(J«)  «  -  XI  *  BTHETA  ♦  X2  *  < THETA!  1%. J%)  ♦  BTHETA) 

390  NG(J%)  ■  ( LM  *  Vt  ♦  .3  •  IP  *  VP)  *  SIM  (THETA< I*. J%) ) 

400  IF  IP  *  COS  (THETA ( !%. J%) )  <  Q  THEN  GOTO  404 

402  MV(J%)  >  .23  •  (UP  •  COS  (THETA  ( 1%.  J%)  >  >  a  2  -  D  a  2)  *  RA  *  CA 
*  DP  *  (UA  a  2) 

403  LS  •  D  /  COS  (THETA( 1%. J%) ) 

404  GOTO  408 
404  MV(J%)  -  0 

407  LS  -  LP 

408  IH  «  IH  ♦  (3.1414  /  12)  *  (OP  A  2)  •  RV  •  CM  •  (AS  a  3) 

410  DCTA  -13/10 

413  SUM  -  0 

420  FOR  K  »  1  TO  18 

423  S  >  DLTA  *  N 

430  X  -  S  *  SIN  <THETA( I*. J%) ) 

433  T  -  -  D  ♦  S  *  COS  (THETA ( 1%. J%) ) 

440  «  -  .3  *  SIGMA  *  KV  *  < <  EXP  ( XA  *  DT  ♦  XA  *  T)  ♦  EXP  (  -  XA  *  DT 

-  XA  *  *))  /  (  EXP  (XA  *  DT)  -  EXP  <  -  XA  *  DT)>)  *  COS  (XA  «  X  -  S 

IGMA  •  T(I%.J%)> 

443  V  •  .3  •  SIGMA  *  MV  •  U  EXP  (XA  •  DT  ♦  XA  *  T)  -  EXP  (  -  XA  •  DT 

-  XA  *  T)>  /  <  EXP  (XA  *  DT)  -  EXP  (  -  XA  *  DT) ) )  *  SIN  (XA  »  X  -  S 

IGMA  «  T(I%.J%)> 

430  UDOT  -  .3  •  (SIGMA  a  2)  *  KV  •  ( (  EXP  (XA  *  DT  ♦  XA  *  T)  ♦  EXP  ( 

-  XA  *  DT  -  XA  *  D)  /  (EXP  (XA  *  DT)  -  EXP  (  -  XA  *  DT) ) )  *  SIN  (X 

A  *  X  -  SIGMA  *  T( II. J%) ) 

433  VDOT  >  -  .3  *  (SIGMA  a  2)  *  HV  «  ((  EXP  (XA  *  DT  ♦  XA  *  T)  -  EXP 

(  -  XA  •  DT  -  XA  *  F>)  /  (  EXP  (XA  •  DT)  -  EXP  (  -  XA  *  DT) ) )  »  COS 

(XA  •  X  -  SIGMA  •  Ttll.JS) ) 

440  FO  -  .3  •  RV  «  CV  •  DP  •  ( <V  *  COS  (THETA ( II .  J%) )  -  V  »  SIN  (THE 

TA(  1%.  J%) )  ♦  VC  »  COS  (THETA ( I%.J%)  >  -  S  *  OMEGA!  !%.<!%))  a  2)  *  SGN 

(U  •  COS  (THETA!  I*.  J*))  -  V  •  SIN  (THETA! I*. J%> >  ♦  UC  •  COS  (THETA( 

I*.J*>>  -  S  •  OMEGA! I*. J*) ) 

442  FI  ■  .23  *  3.1414  •  RV  *  CM  «  (DP  A  2)  •  (UDOT  •  COS  (THETA(I*.J* 
))  -  VDOT  *  SIN.  (THETA(  1%.  J%) )  -  UC  *  (  SIN  (THETA(  IV  J%) ) )  •  OMEGAd 
V J%)> 


OMEGA < 1%. 0 )  +  (H  /  4)  *  <M(0>  ♦  2  *  M(l>  ♦  2  *  M 


THETA  (DEC)  ■  " ; THETA ( 1%. 0)  *  180 


•44  IF  M  >  9.7  THEN  GOTO  444 
444  SUM  »  SUM  *  S  *  <  FD  *  FI) 

447  GOTO  470 

444  SUM  a  SUM  ♦  .5  *  5*  (FD  ♦  FI) 

470  NEXT  M 

440  MC(U%)  -  OITA  *  SUM 

490  M(J%)  a  (  -  MH(J*)  ♦  NC(J%)  ♦  NV(J%)  ♦  MC(J%)>  /  1H 
700  I?  J%  •  3  THEN  GOTO  730 

710  THETA(I%.J%  ♦  1)  a  THETA ( 1%.  0 )  ♦  (T(I».J*  ♦  1)  -  T(I%.0)>  *  X(J%) 
720  OMEGA (  1%.  <J%  *  l)  a  OMEGA! I%.0>  +  <TCI*.J*  ♦  t)  -  T(I«.0>)  •  M(J%) 
730  J%  «  J%  ♦  l 
740  GOTO  310 

730  THETA ( 1%  ♦  1.0)  a  THETA(I%.0)  ♦  <H  /  4)  *  <E(0>  ♦  2  *  X<1>  +  2  *  X 
(2)  +  X ( 3 ) ) 

740  OMEGA ( 1%  ♦  1.0)  a  OMEGA(1«.0)  ♦  (H  /  4)  *  (M(0)  *  1  •  M(l)  ♦  2  *  M 
(2)  ♦  M( 3 ) ) 

770  I*  a  I*  ♦  1 

773  PRINT  “TIME  <S)  a  - :T< 1%. 0) i "  THETA  (DEG)  a  THETA ( 1%. 0)  *  180 

/  3.1414 

780  GOTO  430 

783  PRINT 

790  INPUT  “INPUT  INTERVAL  (•  OF  TIME  STEPS)  FOR  OUTPUT  " ; ITRVL 
800  PR*  1 
410  I*  a  0 

420  IF  THETA ( 1%. 0)  )  (BTHETA)  THEN  GOTO  840 
830  IF  THETA ( 1% . 0 )  ♦  BTHETA  <  0  THEN  GOTO  880 
840  MH( 0 )  a  Xl  i  THETA ( 1% . 0 > 

830  GOTO  890 

440  NH(0)  a  XI  *  BTHETA  ♦  X2  *  <THETAU».0)  -  BTHETA) 

870  GOTO  890 

480  MH(0)  a  -  XI  *  BTHETA  ♦  X2  *  (THETA(I*.0)  ♦  BTHETA) 

890  PRINT  “T  (S)  -  *:T(II.0):*  THETA  (DEG)  -  - ;THETA< I*.0)  »  180  /  3 
.1414;“  HINGE  MOMENT  (FT-LBS)  a  “:MH<0> 

893  PRINT 

900  I*  a  I*  ♦  ITRVL 

910  IF  1%  >  L%  THEN  GOTO  930 

920  GOTO  420 

930  PR*  0 

94  0  END 


COL PILE 


:  PRINT  “PILE  DYNAMICS  DURING  COLLISION” 

S  PRINT 

iO  INPUT  “INPUT  STSPSIZE  (S)  “;H 
IS  PRINT 

20  INPUT  “INPUT  1ST  STIFFNESS  CONSTANT  XI  CFT-L8S/RAD)  “:X1 
23  PRINT 

30  INPUT  “INPUT  2ND  STIFFNESS  CONSTANT  X2  (FT-LBS/RAD)  “:K2 
35  PRINT 

40  INPUT  “INPUT  BREAKPOINT  ANCLE  (DEG)  * ; BA 
45  PRINT 

SO  INPUT  “INPUT  LOAD  WEIGHT  UBS)  “:VL 
33  PRINT 

40  INPUT  “INPUT  PILE  WEIGHT  (LBS)  “;WP 
43  PRINT 

70  INPUT  "INPUT  PILE  LENGTH  <FT>  “;LP 
73  PRINT 

10  INPUT  “INPUT  PILE  DIAMETER  (FT)  “:OP 


as 

PRINT 

70 

INPUT  “INPUT  LENCTH 

TO  LOAD  <FT)  “;LH 

73 

PRINT 

100 

INPUT  ‘INPUT  DEPTH 

TO  HINGE  (FT)  ":D 

105 

PRINT 

110 

INPUT  “INPUT  TOTAL 

DEPTH  (FT)  ";DT 

115 

PRINT 

120 

INPUT  “INPUT  CURRENT  VELOCITY  <FT/3>  “;UC 

122 

PRINT 

125 

INPUT  “INPUT  BARGE 

COEFFICIENT  OF  FRICTION  “;MU 

127 

PRINT 

130 

INPUT  "INPUT  BARGE 

FREEBOARD  (FT)  “ ; FB 

133 

PRINT 

140 

INPUT  “INPUT  BARGE 

DRAFT  (FT)  “ ;DS 

143 

PRINT 

130 

INPUT  "INPUT  BARGE 

BOV  ANGLE  (DEG)  " ;FA 

155 

PRINT 

140 

INPUT  "INPUT  BARGE 

LENCTH  (FT)  ":L8 

143 

PRINT 

170 

INPUT  “INPUT  BARGE 

SPEED  (FT/S)  " ; UB 

173 

PRINT 

ISO 

CV  a  1.0 

113 

RV  a  44.0  /  32.2 

170 

CM  a  2.0 

173 

PI  a  3.1414 

200 

BTHETA  a  BA  *  PI  / 

110 

210 

FTHETA  a  FA  *  PI  / 

110 

220  LF  «  <08  ♦  FS>  *  TAN  (FTHETA) 

223  PHI  a  ATN  (MU) 

230  LX  >  LB  -  LF 
240  Hi  •  0  •  SI 

230  TF  .  <<0  ♦  FB)  /  UB)  •  TAN  (FTHETA) 
240  TNXT  a  SQR  (IIP  a  2  /  KB  a  2)  .  1) 
270  KTHXTA  a  ATN  (TNXT! 

210  TX  -  CHB  *  TNXT  ♦  LF)  /  UB 
2P0  TX  •  TX  ♦  LX  /  UB 
300  ML  •  WL  /  32.2 
310  MP  ■  WP  /-  32.2 
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320  IM  a  (MP  *  LP  A  2)  /  3  +  ML  *  LM  a  2 
330  T  «  S 

333  PR#  1 

34  0  PRINT  "IMPACT  AT  TOP  OF  3 A AGE  BOV" 

330  PRINT  “TIME  (S)  a  0  THETA  (DEC)  -  0  " 

355  IH  a  IM  ♦  RV  *  CM  *  PI  *  (DP  a  2)  ■  (D  a  ])  M2 
340  MDT  ■  IH  *  UB  /  (Of  FB) 

370  FD7  «  (IH  *  UH>  /  ((D  +  FS>  a  2> 

375  MA  a  CM  *  RV  «  PI  *  (DP  a  1)  *  D  /  4 

30 0  HOT  a  (IH  /  (D  ♦  FB)  -  (ML  *  LM  ♦  .5  *  (MP  «  LP  +  MA  *  D>>>  *  VB  / 
(D  +  FB) 

390  VOT  a  0 

400  PRINT  "BARGE  MOMENT  IMPULSE  (FT-LB-S)  a  " ; MDT ; "  BARGE  FORCE  IMPU 
L3E  (LB-S)  a  " ; FDT 

410  PRINT  “H0R2  REACT  IMPULSE  (LB-S)  a  " ; HDT ; "  VERT  REACT  IMPULSE  (L 
B-S )  a  " ;VDT 
415  PRINT 

410  IF  (T  ♦  H)  )  »  TF  THEN  GOTO  390 

420  T  a  T  ♦  H 

430  THETA  a  ATN  (UB  *  T  /  (D  ♦  FB)) 

431  IF  (D  ♦  FB)  /  COS  (THETA)  (  a  LP  THEN  GOTO  440 

432  THETA  a  ATN  ( ( 1  -  (D  +  FB)  /  (LP  *  COS  (THETA)))  *  TAN  (FTHETA) 
♦  UB  *  T  /  (LP  *  COS  (THETA))) 

433  OMECA  a  UB  /  (L?  *  (  COS  (THETA)  ♦  TAN  (FTHETA)  *  SIN  (THETA))) 

434  ALPHA  a  (  SIN  (THETA)  -  TAN  (FTHETA)  *  COS  (THETA))  *  ((UB  I  LP) 
a  2)  /  ((  COS  (THETA)  ♦  TAN  (FTHETA)  *  SIN  (THETA))  a  3) 

435  LC  a  LP 

434  GOTO  440 

440  OMEGA  a  UB  *  ( (  COS  (THETA))  a  2)  /  (D  ♦  FB) 

430  ALPHA  a  -  2  *  ((UB  /  <D  ♦  FB))  a  2)  *  ((  COS  (THETA))  a  3)  •  SIN 

(THETA) 

435  LC  a  (D  ♦  FB)  /  COS  (THETA) 

440  PRINT  "TIME  (S)  .  «:T;“  THETA  (DEG)  -  THETA  *  ISO  /  PI 

440  IF  LP  *  COS  (THETA)  (  D  THEN  GOTO  474 

470  LS  a  D  /  COS  (THETA) 

472  GOTO  400 
474  LS  a  LP 

400  IH  a  IM  ♦  RV  •  CM  *  PI  *  (DP  A  2)  *  (LS  A  3)  /  12 

<05  MA  a  CM  *  RV  t  pi  *  (DP  a  2)  *  LS  /  4 

4)0  COSUB  2000 
300  GOSUB  3000 

510  H  a  IH  •  ALPHA  ♦  MH  -  LM  «  VI  *  SIN  (THETA)  -  .5  *  LP  »  VP  *  SIN 
(THETA)  -  MC 
SIS  COSUB  4000 

520  IF  (D  ♦  FB)  /  COS  (THETA)  >  LP  THEN  GOTO  552 

330  P  a  M  /  (LC  •  COS  (PHI)) 

540  RH  a  F  •  COS  (THETA  -  PHI )  *  FC  *  COS  (THETA)  ♦  (ML  *  LM  ♦  . 5  * 

MP  *  LP)  *  ((OMEGA  a  2)  *  SIN  (THETA)  -  ALPHA  *  COS  (THETA))  -  MA  * 

.3  *  IS  *  ALPHA  •  COS  (THETA) 

350  RV  a  F  *  SIN  (THETA  -  PHI)  ♦  FC  *  SIN  (THETA)  ♦  VL  ♦  VP  -  (ML  * 

LM  ♦  .5  *  MP  *  LP)  *  (ALPHA  «  SIN  (THETA)  ♦  (OMEGA  a  2)  *  COS  (THETA 

))  -  MA  •  .3  *  LS  *  ALPHA  *  SIN  (THETA) 

331  GOTO  340 

332  F  a  M  /  (LP  *  COS  (THETA  -  FTHETA  ♦  PHD) 

533  RH  a  F  t  COS  (FTHETA  -  PHI)  ♦  FC  *  COS  (THETA)  ♦  (ML  *  LM  ♦  .5  * 
MP  •  LP)  «  ((OMEGA  a  2)  •  SIN  (THETA)  -  ALPHA  «  COS  (THETA))  -  MA  * 
.3  •  IS  *  ALPHA  •  COS  (THETA) 

334  RV  a  F  *  SIN  (FTHETA  -  PHI )  ♦  FC  *  SIN  (THETA)  t  VL  +  VF  -  (ML  * 

LM  ♦  .3  •  MP  *  LP)  *  (ALPHA  »  SIN  (THETA)  ♦  (OMEGA  a  2)  *  COS  (THET 

A) )  -  MA  «  .3  *  LS  *  ALPHA  *  SIN  (THETA) 
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540  PRINT  "BARGE  MOMENT  <?T-tBS>  -  ":M;"  BARGE  FORCE  (LBS)  *  ":F 

570  PRINT  “HINGE  MOMENT  (FT-LBS)  «  ":MH:“  HORZ  REACT  (LBS)  a  “ ; RH : " 

VERT  REACT  (LBS)  a  “;RV 
575  PRINT 
580  GOTO  411 

59 0  PRINT  “IMPACT  AT  BOTTOM  OP  BARGE  BOW" 

400  PRINT  "TIME  (S>  a  “;TF;"  THETA  (DEC)  a  " ; PTHETA  *  180  /  PI 

402  IP  LP  *  COS  (PTHETA)  <  O  THEN  GOTO  405 

403  LS  a  O  /  COS  (PTHETA) 

404  GOTO  404 

405  LS  •  LP 

404  IH  «  IM  ♦  RW  •  CM  *  PI  «  (DP  a  2)  «  (LS  a  3)  /  12 

408  MA  a  CM  *  RV  *  PI  *  (DP  a  I)  •  LS  /  4 

410  LC  -  HB  /  COS  (PTHETA) 

420  MDT  a  IH  *  UB  *  ((  COS  (PTHETA))  a  2)  *  (1  /  HB  -  < 1  /  (D  +  PI) ) ) 
430  POT  a  MDT  /  LC 

44  0  ROT  a  (IH  /  LC  -  (ML  *  LM  ♦  .  5  *  (MP  *  LP  ♦  MA  *  LS ) ) )  *  VB  *  ( (  C 
OS  (PTHETA))  a  2)  •  ((l  /  HB)  -  (1  /  (D  ♦  F3> ) ) 

450  HOT  a  RDT  *  COS  (PTHETA) 

440  VDT  a  RDT  *  SIN  (PTHETA) 

470  PRINT  "BARGE  MOMENT  IMPULSE  (FT-LS-3)  a  »:MDT:»  BARGE  FORCE  IMP 
ULSE  (LB-S)  a  " ; FDT 

430  PRINT  "HORZ  REACT  IMPULSE  (LB-S)  a  " ;HDT: "VERT  REACT  IMPULSE  (LB- 
S)  a  " ; VDT 
485  PRINT 

488  IP  (T  ♦  H)  >  a  TX  THEN  GOTO  710 
700  T  -  T  ♦  H 

710  THETA  a  ATN  ( (UB  «  T  -  LP)  /  HB) 

720  OMEGA  .  (UB  /  HB)  «  ((  COS  (THETA))  a  2) 

730  ALPHA  a  -  2  *  ((UB  /  HB)  a  2)  *  ((  COS  (THETA))  a  3)  *  SIN  (THET 

A) 

740  LC  a  HB  /  COS  (THETA) 

750  IF  LP  *  COS  (THETA)  <  D  THEN  GOTO  780 
740  LS  a  D  /  COS  (THETA) 

770  GOTO  770 
780  LS  a  L? 

770  IH  a  IM  ♦  RW  *  CM  *  PI  «  (DP  a  2)  *  (LS  A  3)  f  12 

795  MA  a  CM  *  RW  *  PI  *  (DP  a  2)  *  LS  /  4 

800  PRINT  "TIME  (S)  a  ";T;"  THETA  (DEG)  a  " ; THETA  *  180  /  PI 
810  COSUB  2000 
120  COSUB  3000 

830  M  a  IH  •  ALPHA  ♦  MH  -  LM  •  WL  *  SIN  (THETA)  -  .5  *  LP  *  WP  *  SIN 
(THETA)  -  MC 

140  F  a  M  /  (LC  •  COS  (PHD) 

850  GOSUB  4000 

840  RH  a  p  *  COS  (THETA  -  PHI)  ♦  PC  *  COS  (THETA)  ♦  (ML  *  LM  ♦  .5  * 

MP  *  LP)  *  ((ONECA  a  2)  «  SIN  (THETA)  -  ALPHA  *  COS  (THETA))  -  MA  * 

.5  *  LS  •  ALPHA  *  COS  (THETA) 

870  RV  a  p  »  SIN  (THETA  -  PHI)  ♦  PC  *  SIN  (THETA)  ♦  WL  ♦  WP  -  (ML  * 
LM  ♦  .5  *  MP  *  LP)  *  (ALPHA  «  SIN  (THETA)  ♦  (OMEGA  a  2)  *  COS  (THETA 
>)  -  HA  *  .5  •  LS  *  ALPHA  *  SIN  (THETA) 

880  PRINT  "BARGE  MOMENT  (PT-LBS)a  ";H;"  BARGE  FORCE  (LlS)a  «  ;F 

870  PRINT  "HINGE  MOMENT  (FT-LBS) a  " ;NH ; "  HORZ  REACT  (LBS) a  ";RH;”  V 

ERT  REACT  (LBS)  a  '•  .  RV 

875  PRINT 

700  GOTO  438 

710  PRINT  "PILE  IN  CONTACT  WITH  BARGE  BOTTOM" 

720  PRINT  "TIME  DURATION  (S>  IS  FROM  "iTX;"  TO  " ;TR; 

725  PRINT  "THETA  (DEC)  a  " ; XTHETA  *  180  /  PI 
730  THETA  a  XTHETA 
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440  OMEGA  *  0 
930  ALPHA  •  0 
9*0  LC  •  L? 

970  LS  •  L? 

900  COSUB  2000 
990  COSUB  0000 

1000  M  a  MH  -  LM  «  VL  *  SIN  < THETA)  -  . S  *  IP  *  VP  *  SIN  (THETA)  -  M 
C 

1010  F  a  M  /  <LP  *  SIN  (XTHETA  ♦  PHD) 

1030  COSUB  4000 

1030  RH  «  F  *  SIN  (PHI)  ♦  PC  *  COS  (THETA) 

1040  IV  a  M  COS  (PHI)  ♦  FC  *  SIN  (THETA)  ♦  VI  *  VP 

10S0  PRINT  "BARGS  MOMENT  (FT-LBS)  a  «;M;"  BARGE  FORCE  (LBS) ■  ":F 

10*0  PRINT  "HINCE  MOMENT  (FT-LBS)  a  " ; MH ; "  HORZ  REACT  (LBS)  a  " ; RH ; “ 

VERT  REACT  (LBS)  a  ":RV 

1070  PRINT  "PILE  RECOVERY  IS  INITIATED  AT  T  <S>  a  “ ;TR 
1000  PRO  0 
1090  END 

3000  IF  THETA  >  ( BTHETA )  THEN  GOTO  3040 
3010  IF  THETA  (  (  -  BTHETA)  THEN  COTO  30*0 
3020  MH  a  XI  «  THETA 
3030  COTO  2070 

3040  MH  a  XI  *  BTHETA  ♦  X2  *  (THETA  -  BTHETA) 

30S0  COTO  2070 

20*0  MH  -  -  XI  «  BTHETA  ♦  X3  •  (THETA  ♦  BTHETA) 

3070  RETURN 

3000  0LTA  a  LS  /  10 

3010  SUM  a  0 

3020  FOR  N  a  1  TO  10 

3030  S  a  DLTA  *  N 

3040  FD  a  .5  «  RW  *  CV  •  DP  *  ( (UC  •  COS  (THETA)  -  S  *  OMEGA)  a  2)  • 
3CN  (UC  *  COS  (THETA)  -  S  *  OMEGA) 

3030  FI  a  -  .25  *  PI  *  RW  *  CM  *  (DP  a  2)  *  UC  *  OMEGA  *  SIN  (THETA) 

30*0  IF  N  >  9.9  THEN  GOTO  3090 
3070  SUM  a  SUM  ♦  S  *  (FD  a  FI) 

3000  COTO  3100 

3090  SUM  a  SUM  ♦  .3  *  S  *  (FD  ♦  FI) 

3100  NEXT  N 

3110  MC  a  DLTA  *  SUM 

3120  RETURN 

4000  DLTA  a  LS  /  10 

4010  SUM  a  0 

4020  FOR  N  a  0  TO  10 

4030  S  a  DLTA  •  N 

4040  FD  a  .5  •  RW  «  CV  «  OP  «  ((UC  *  COS  (THETA)  -  S  •  OMEGA)  a  2)  • 
3GN  (UC  •  COS  (THETA)  -  S  *  OMEGA) 

4050  FI  a  -  .25  *  PI  *  RW  *  CM  *  (DP  a  2)  *  UC  *  OMECA  *  SIN  (THETA) 

40*0  IF  N  (  .9  THEN  GOTO  4100 
4070  IF  N  )  9.9  THEN  GOTO  4100 
4000  SUM  -  SUM  ♦  .  S  *  (FD  ♦  FI) 

4090  COTO  4110 

4100  SUM  a  SUM  ♦  .3  *  (FD  ♦  FI) 

4110  NEXT  N 

4120  FC  a  DLTA  *  SUM 

4130  RETURN 
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RECPILE 


A 

PRINT  “ 

PILE  MOTION  DURING  RECOVERY" 

5 

PRINT 

10 

INPUT 

"INPUT  STEPS I ZE  (S>  ";H 

15 

PRINT 

:o 

INPUT 

"INPUT  MAX  T  <S>  " ;MAX 

23 

PRINT 

30 

L%  a  INT  (MAX  /  H) 

| 

40 

INPUT 

“INPUT  1ST  ST1PFNE5S  CONSTANT 

XI  (FT-LBS/RAD) 

" ;  XI 

05 

PRINT 

SO 

INPUT 

"INPUT  2ND  STIFFNESS  CONSTANT 

K2  (FT- LBS /RAD) 

“:K2 

vv; 

S3 

PRINT 

•  0 

INPUT 

"INPUT  BREAKPOINT  ANGLE  (DEG) 

" :  BA 

01 

BTKETA 

a  BA  *  3.1410  /  180 

i 

05 

PRINT 

70 

INPUT 

"INPUT  LOAD  WGT  (LBS)  “ ;WL 

?  ‘ 

73 

PRINT 

10 

ML  •  VL  /  32.2 

70 

INPUT 

“INPUT  PILE  WGT  (LBS)  “;VP 

m  *v 

75 

PRINT 

1 

100  HP  a  ' 

UP  /  32.2 

110  INPUT  “INPUT  PILE  LENGTH  <PT)  " iLP 
IIS  PRINT 

110  INPUT  "INPUT  PILE  DIAMETER  (FT)  " ;DP 
123  PRINT 

130  INPUT  "INPUT  LENCTH  TO  LOAD  (PT)  "  jIM 
133  PRINT 

130  INPUT  "INPUT  DEPTH  TO  HINGE  (PT)  " ;D 
153  PRINT 

UO  INPUT  "INPUT  TOTAL  DEPTH  (PT)  " ;DT 
US  PRINT 


100 

INPUT  "INPUT  CUR  VEL  (FT/S)  ";UC 

IIS 

PRINT 

170 

INPUT  "INPUT  INITIAL  THETA  (DEG)  “;T0 

173 

PRINT 

200 

XTHSTA  a  TO  *  3.1410  /  180 

230 

RV  a  04.0  /  32.2 

210 

CV  a  1.0 

1 

300 

CM  a  2.0 

310 

IM  a  (HP  •  LP  A  2)  /  3  ♦  ML  *  LM  A  2 

320 

DIM  T(L%  ♦  2.3) 

330 

DIM  THETA ( L%  ♦  2.3) 

340 

DIM  OMEGA ( L*  ♦  2.3) 

K-[- 

330 

DIM  X( 3 ) 

300 

DIM  M( 3 ) 

*  *  „ 

370 

DIM  MHO) 

380 

DIM  MG( 3) 

•*  *- 

400 

DIM  MC( 3 ) 

410 

I*  ■  0 

i 

420 

T(0.0>  a  0 

430  THETA (0.0)  a  XTKCTA 

007  PRINT  "TINE  (S>  a  ";T(I*.0>;"  THETA  (DEG)  a  " ;THETA( I*. 0)  «  110 
/  3.1010 

030  T(I*  ♦  1.0)  a  T(I*.0>  ♦  H 

000  If  T(I%  ♦  1.0)  >  MAE  THEN  GOTO  705 

070  TCl.l)  a  T(1«.0)  ♦  .3  *  H 
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.•  .  * 
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483  Tv I%.2)  »  T( I*. i; 

490  7(1%. 3)  ■  T( :%. 0)  ♦  H 
S00  J%  a  0 

310  XCJ%>  a  OMEGA ( 1%. J%> 

320  IP  THETAC 1%. J%>  >  C BTKETA)  THEN 
330  If  THETA < 1%. J%)  *  3 THETA  C  0  TH: 
340  MH(J%)  -  XI  *  THETA ( 1%. J%) 

330  GOTO  390 

340  MH( J%>  .  XI  «  BTHETA  a  X2  *  (THE' 
370  GOTO  390 

380  MH(J%)  a  -  XI  *  STHSTA  a  X2  *  C‘ 
390  MGC J%)  a  ( LJ1  *  VL  a  .5  *  IP  *  VP 


IEN  GOTO  S  40 
THEN  GOTO  380 


< THETA < !%.J%>  -  BTHETA) 


•  ( THETA < 1%. J%)  a  BTHETA) 

•  VP)  *  SIN  (THETAC 1%. J%) > 


400  IF  LP  * 

403  LS  a  O  / 

404  GOTO  408 
404  LS  a  LP 
408  IH  •  IM  ♦ 
410  OLTA  a  LS 
413  SVM  a  a 
420  FOR  N  a 
423  S  a  0LTA  < 
440  F0  a  .3  • 
,J%>)  a  2)  * 


COS  <TH£TA< 1%. J%) )  <  0  THEN 
COS  (THETAC !%. J%) > 


(3.  1414 
/  10 

TO  10 
N 

RV  *  CV 
3CN  CUC 


GOTO  404 


442  FI  a  -  .23  * 
IN  (THETAC 1%. J%) ) 


3.1414 


>F  *  (CUC  •  COS  CTHETAC I%.J%> >  -  S  *  OMEGA! 1% 
COS  (THETAC 1%.J%> )  -  S  *  OMECAC 1%. J%> ) 

«  RV  *  CM  *  COP  a  2)  *  UC  *  OMEGA C 1%. J%)  »  S 


MC  C  J% ) ) 


CTC 1% , J% 
CTC 1%.  J% 


>  TC 1%. 8)  ) 
-  TC 1%. 0 > ) 


444  IF  N  )  9.9  THEN  GOTO  441 
444  SUM  a  SUM  a  S  *  CFO  ♦  FI) 

447  GOTO  470 

448  SUM  a  SUM  a  .3  *  S  *  CFO  a  FI) 

470  NEXT  N 

480  MCC J%)  a  OLTA  *  SUM 

490  M(J%)  a  C  -  MHCJ%)  a  MGCJ%>  *  MCCJ%)>  /  IH 
700  IF  J%  a  3  THEN  GOTO  730 

710  THETAC 1%.  J%  a  1)  a  THETACI%.0>  ♦  CTCI%.J%  a  1)  -  T(I%.0>)  • 
720  OMECAC 1%. J%  a  1)  .  OMEGACI%.0)  a  CT(I%.J%  a  l>  -  TCI%.0>)  * 
730  J%  a  J%  ♦  1 
740  GOTO  SlO 

750  THETAC 1%  ♦  1.0)  a  THETAC1%.0)  a  CK  /  4)  *  CXCO)  *  2  •  XC1)  < 
C  2 )  a  KC3) ) 

740  OMECAC 1%  a  1.0)  a  OMEGACI%.0)  a  CH  /  4)  *  CMCO)  a  2  •  MCI)  * 
C  2  >  a  MC  3 ) ) 

770  1%  a  :%  a  1 

773  PRINT  -TIME  CS)  a  -;TCI%.0);-  THETA  COEG)  a  " -.THETAC  1%.  0) 

/  3.1414 
780  GOTO  430 
783  PRINT 

790  INPUT  -INPUT  INTERVAL  CO  OF  TIME  STEPS)  TOR  OUTPUT  " ; ITRVL 
800  PRO  1 


XC  J%> 
MCJ%> 


a  THETAC I%.0) 


OMEGA C 1%. 0 ) 


CXCO) 


(MC  0) 


*  MCI) 


a  -  ;T(  1%.  0)  ;  "  THETA  COEG)  a  - -.THETAC  1%.  0) 


/  3.  14 
78  0  C< 

783  P 
790  II 
800  P 
810  1% 
120  I 
830  I 


IP  THETAC I%.0> 
IF  THETAC I%,0) 


>  (BTHETA)  THEN  GOTO  840 
a  BTKETA  C  0  THEN  GOTO  880 


840  MHCO)  a  XI  •  THETAC 1%. 0 ) 

830  GOTO  890 

140  MHCO)  a  XI  •  BTHETA  a  X2  « 
870  GOTO  890 

880  MH<  0)  a  -  XI  «  BTHETA  a  X2 
890  PRINT  -T  CS)  a  -;TCI%.0);" 
.1414;-  HINGE  MOMENT  (FT-LBS) 
893  PRINT 
900  1%  a  1%  a  ITRVL 


CTHETAC 11.0)  -  BTHETA) 


•  CTHETAC I%.0> 
THETA  (DEG)  a 
- ;MH( 0) 


a  BTHETA) 

" ; THETAC 1%. 0 ) 


100  / 


XIII.  APPENOIX  B:  OTHER  HINGE  CONCEPTS 


The  rubber  tube  analysis  results  indicated  that  it  would  be  necessary  to 
consider  hinge  components  which  make  use  of  springs.  Many  spring  loaded 
mechanisms  were  investigated  but  not  adopted.  Several  of  the  concepts  which 
warranted  serious  consideration  are  presented  in  this  appendix.  In  each  case, 
the  reason  for  their  ultimate  rejection  is  discussed. 

A  very  simple  concept  using  a  central  spring  and  a  circular  base  is  shown 
in  Fig.  31.  The  spring  is  housed  in  the  hollow  pile,  and  the  lower  part  of  the 
pile  is  widened  for  increased  stability.  When  tipped,  the  pile  rocks  on  the  base 
contact  point,  and  spring  tension  provides  the  restoring  moment.  The  spring  is 
prestressed  so  the  restoring  moment  develops  at  very  small  angle  changes.  While 
promising,  this  particular  concept  was  rejected  because  there  is  no  simple  way 
to  ensure  that  base  contact  Is  maintained  during  lower  bow  rake  collision  impact. 

Peripheral  stay/point  universal  joint  concepts  were  then  developed  to 
provide  a  fixed  point  of  attachment  as  well  as  obtain  Initial  stiffness  by  means 
of  prestressing.  As  shown  In  Fig.  32,  the  pile  Is  centrally  supported  using  a 
crossed  axis  universal  joint,  and  its  angular  motion  is  restrained  by  stays 
placed  on  the  periphery.  The  arrangement  is  such  that  If  there  Is  an  angular 
change  from  the  vertical,  the  stay  on  the  outside  of  the  bend  immediately  sustains 
the  entire  prestressing  force.  Thus  a  large  restoring  moment  Is  generated  at 
very  small  angles.  This  concept  was  not  used  In  the  manner  shown  in  Fig.  32 
because  the  joint  requires  90  deg  of  rotation  which  Is  beyond  the  range  of  the 
standard,  commercially  available  U-joint.  In  addition,  the  stay  moment  arm  with 
respect  to  the  joint  decreases  at  large  angles.  Consequently  the  restoring 
moment  is  reduced  below  that  required  to  return  the  pile  to  the  upright  position 
from  a  full  knockdown. 

The  segmented  pile/internal  stay  hinge  system  shown  In  Fig.  33  provides  full 

articulation.  Stay  moment  arm  with  respect  to  the  hinge  axes,  however,  Is  further 

B-l 


Fig.  31,  Central  spring/circular  base  hinge  concept.  F 
is  the  force  due  to  a  central,  Internal  "sprlnf 
which  Is  attached  to  the  base  (circular  in  plan 
view)  by  the  cable  shown. 


Fig.  32.  The  peripheral  stay' /point  universal  joint  concept. 

Four  stays  are  used  which  are  prestressed  by  F  according 
to  Eq.  12.  s 


Fig.  33.  The  segmented  plle/internal  stay  concept.  Four 

Internal  stays  are  guided  to  lead  across  the  joint 
notches.  The  stays  are  prestressed  by  F  . 


reduced.  The  concept  was,  therefore,  rejected  due  to  inadequate  restoring 
moment  stiffness. 

The  last  concept  presented  here  involves  the  use  of  a  break- away  shear  pin. 
The  spring  in  any  of  the  previously  described  systems  is  replaced  by  a  tension 
member.  This  member  is  prestressed  and  pinned  at  the  top  as  shown  in  Fig.  34. 
When  the  pile  is  struck,  the  pin  will  fail  in  shear  and  the  pile  rotates  freely 
at  the  hinge.  The  pile  system  remains  connected  enabling  its  position  to  be 
easily  relocated.  The  system  can  then  be  reassembled  without  the  use  of  a  pile 
driver  equipped  barge.  Though  simple  and  workable,  this  system  does  not  provide 
an  automatic  return  to  the  operating  position.  Because  assistance  is  required, 
this  concept  was  judged  not  able  to  meet  the  fundamental  design  objective. 

While  the  hinge  concepts  described  in  this  appendix  were  not  used,  they 
were  not  without  merit.  In  fact  the  final  hinge  concept  adopted,  as  shown  in 
Fig.  5,  evolved  from  these  earlier  efforts  and  incorporates  several  of  their 
more  useful  features. 


XIV.  APPENDIX  C:  SPRING  OPTIONS 


Two  principal  concepts  for  creating  the  "spring"  force  F$  are  possible  - 
an  airbag/cyl inder  concept  and  an  elastic  member  such  as  a  rubber  band  or  an 
array  of  nylon  ropes.  It  was  found  that  helical  steel  springs  are  not  suitable. 
Helical  springs  which  possess  the  required  stiffness  and  elastic  range  are  much 
too  large  to  fit  inside  the  pile. 

In  the  airbag  concept,  shown  in  Fig.  35,  the  bag  is  contained  between  a 
platform  fixed  to  the  inside  of  the  pile  wall  and  a  moveable  piston.  The  piston 
rod  lies  along  the  pile  centerline  and  passes  through  the  platform  to  a  point 
where  the  stays  are  attached.  The  airbag  is  filled  from  the  surface  to  a  pressure 
corresponding  to  the  necessary  prestress  force  Fs>  When  the  pile  moves  from  the 
vertical,  the  air  is  further  compressed  contributing  to  a  positive  k2  value. 

The  bag  operating  pressure  is  in  the  vicinity  of  400  psi  which  allows  easy  in¬ 
flation.  Conversations  with  engineers  familiar  with  commercially  available  air¬ 
bags,  however,  suggest  that  this  arrangement  pushes  the  limit  of  current  technology. 

The  other  acceptable  spring  concept  involves  the  use  of  an  elastic  member. 

The  stays  are  connected  to  one  end  of  the  member,  and  the  member  is  then  stretched 
until  the  load  equals  F$.  A  worm  gear  jack  arrangement,  as  illustrated  in  Fig. 

36,  could  be  used  to  accomplish  the  prestressing.  The  load-deflection  character¬ 
istics  of  the  member  will  be  measured  before  installation  so  that  the  amount  of 
elongation  required  to  produce  the  desired  Fg  is  known.  The  length  increase  at 
F$  should  not  be  the  maximum  elongation  so  the  member  can  stretch  further  during 
collision  without  exceeding  the  elastic  range. 
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Fig.  35.  The  airbag  "spring"  concept.  The  doughnut 

shaped  bag  is  used  to  seal  the  piston-cylinder 
system  shown. 
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Fig.  36.  Worm  gear  jack  for  prestressing  an  elastic 
spring  element. 
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